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Abstract 

One of the main goals in quantum circuit optimisation is to reduce the number of ancillary 
qubits and the depth of computation, to obtain robust computation. However, most of known 
techniques, based on local rewriting rules, for parallelising quantum circuits will require the 
addition of ancilla qubits, leading to an undesired space-time tradeoff. Recently several novel 
approaches based on measurement-based quantum computation (MBQC) techniques attempted 
to resolve this problem. The key element is to explore the global structure of a given circuit, 
defined via translation into a corresponding MBQC pattern. It is known that the parallel power 
of MBQC is superior to the quantum circuit model, and hence in these approaches one could 
apply the MBQC depth optimisation techniques to achieve a lower depth. However, currently, 
once the obtained parallel pattern is translated back to a quantum circuit, one should either 
increase the depth or add ancilla qubits. In this paper we characterise those computations where 
both optimisation could be achieved together. In doing so we present a new connection between 
two MBQC depth optimisation procedures, known as the maximally delayed generalised flow 
and signal shifting. This structural link will allow us to apply an MBQC qubit optimisation 
procedure known as compactification to a large class of pattern including all those obtained from 
any arbitrary quantum circuit. We also present a more efficient algorithm (compared to the 
existing one) for finding the maximally delayed generalised fiow for graph states with fiow. 



*The first two authors have contributed equally to this work. 
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1 Introduction and background material 



We are slowly reaching the classical limit of current technology in decreasing the size of computer 
chips. Hence to avoid quantum effect and also to reduce heat generation parallel computing has 
become the dominant paradigm in computer architecture. On the other hand in the quantum domain 
reducing the depth would be essential for keeping the computation coherent to avoid classical effect. 
In both scenarios designing parallel circuit remains a challenging task. Hence one attempts to 
rewrite locally a given circuit to reduce the depth of the computation. Usually such an approach 
will add ancilla registers to achieve the parallelisation, this being an undesired effect in the quantum 
setting as addition of ancilla qubits could increase the decoherence, breaching the initial purpose 
of parallelisaiton. In this paper we present a general optimisation techniques for quantum circuit 
that exploits the global structure of a given computation to achieve parallelisation with no ancilla 
addition. 

Our technique is based on the translation of a given quantum circuit [1, 2J into a measurement- 
based quantum computation (MBQC) [31(4]. These two models utilise remarkably different informa- 
tion processing tools: while the former is based on unitary evolution of an initially non-entangled 
set of qubits, the latter needs an initial highly- entangled multi-qubit state, where the information 
processing is driven by measurements only. Naturally since the two aforementioned models use 
different information processing tools, hence each model has its own optimisation techniques. In 
the MBQC model, for instance, most of the optimization techniques are based on the identification 
of a more efficient correction structure that is directly linked to the geometry of the underlying 
global entanglement structure. Examples of these techniques are the signal shifting [4J and the 
generalised flow [5j, both discussed in following sections. The so-called standardisation procedure 
[4] can also reduce the number of computational steps by rearranging the MBQC operations into a 
normal form. Moreover, all Pauli measurements in this model can be performed in the beginning of 
the computation [6J, which is a surprising difference from the quantum circuit model. 

On the other hand, most optimization techniques for quantum circuits are based on template 
identification and substitution. For instance in [7J, some circuit identities are used to modify the 
teleportation and dense coding protocols, with the purpose of giving a more intuitive understanding 
of those protocols. Similarly in [8J and [9J a set of circuit identities for reducing the number of gates 
in the circuit for size optimization was given. In contrast to that, in [lOj a useful set of techniques 
for circuit parallelisation was provided, where the number of computational steps is reduced by using 
additional resources. However, as noted in [8J, all the aforementioned circuit optimization techniques 
are basically exchanging a sequence of gates for a different one without any consideration on the 
structure of the complete circuit being optimised. The translation into MBQC would allow us to 
explore the global structure of a given circuit. 

The first such a scheme by back and forth translation between the two models was presented in 
However the backward translation into the circuit required the addition of many ancilla qubits. 
On the other hand several recent works presented an optimised translation scheme from MBQC into 
the circuit where all the non-input qubits are removed, referred hereinafter as compact translation 
[H 1121 [El EH]. The price for doing that happened to be the loss of the optimal depth of the original 
MBQC. The key result of our paper is a new scheme for the parallelisation where the obtained 
MBQC pattern could be translated compactly. Our scheme is based on a new theoretical connection 
between two MBQC depth optimisation procedures, known as the maximally delayed generalised 
flow \T4\ and signal shifting^ which could lead to other interesting observation about MBQC, beyond 
the purpose of this paper as we discuss later. Also our result highlights the fundamental role of the 
MBQC Pauli optimisation in obtaining MBQC parallel structure beyond anything obtainable in the 
quantum circuit model, described later. We prove how our proposed scheme is more optimal in both 
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depth and space compare to the scheme in [TT]. We conclude with a new algorithm for finding the 
maximally delayed generalised flow for graphs with flow with O(n^) steps compared to the exciting 
algorithm in [14j, where n is the number of the nodes in the graph. 

1.1 The MBQC Model 

We review the basic ideas behind the measurement-based quantum computation, with special 
attention to its description in terms of the formal language known as Measurement Calculus [4], and 
the flow theorems [151 IS] • 

In 1999 Chuang and Gottesman described how one could apply arbitrary quantum gates using 
an adaptation of the quantum teleportation model |T6]. This approach was further developed by 
other researchers [TTl IHl [IHl EO] , enabling one in principle to perform arbitrary computations given 
a few primitives: preparation of maximally entangled systems of fixed, small dimension; multi-qubit 
measurements on arbitrary set of qubits; and the possibility of adapting the measurement bases 
depending on earlier measurement outcomes. 

These models of computation draw on measurements to implement the dynamics, and as such 
named collectively the measurement-based model of quantum computation (MBQC), for an overview 
see the paper by Jozsa [21j. An MBQC model using only single qubit measurements was proposed by 
Raussendorf and Briegel in 2001, which became known as the one-way model [6J. The one-way model 
achieves universality through the preparation of a special type of entangled states, the so-called 
cluster states [22j. These states are created with the CZ gate acting on qubits prepared in the state 
1+) = ~^ 1-^)) arranged in a regular lattice, usually the two-dimensional ones. This can be 

relaxed to create more general states with the same interaction over general graphs, creating the 
so-called graph states [23j. Both cluster and graph states can be represented graphically, using 
vertices denoting the qubits and edges for the two-qubit entangling gate CZ. Therefore, the entangled 
resource for the one-way model can be fully represented as graphs. Although two-dimensional cluster 
states can be used as resource for universal quantum computation in the one-way model, arbitrary 
graph states may, or may not, serve for the same purpose; investigating which kinds of entangled 
states are useful resources for MBQC is an active area of research [24l [25l l26t [27]. 

A formal language to describe in a compact way the operations needed for the one-way model 
was proposed in [4J. The language could be easily adapted to any other type of measurement-based 
model hence in the rest of this paper we refer to the general MBQC term instead of the specific 
one-way model as our scheme could be applicable to any MBQC models. In this framework every 
MBQC algorithm (usually referred to as an MBQC pattern) involves a sequence of operations such 
as entangling gates, measurements and feed- forwarding of outcome results to determine further 
measurement bases. A measurement pattern^ or simply a pattern, is defined by a choice of a set of 
working qubits V ^ a subset of input qubits (/), another subset of output qubits (O), and a finite 
sequence of commands acting on qubits in V . Therefore, we consider patterns associated to the 
so-called open graphs: 

Definition 1 (open graph). An open graph is a triplet {G,I,0), where G — iV^E) is a undirected 
graph, and I,0 CV are respectively called input and output vertices. 

An example of an open graph is shown in Figure [3| There are four types of commands, the first 
is the qubit initialisation command Ni that prepares qubit i in the state |+). The input qubits are 
already given as a prepared state. The entangling command Eij = CZij corresponds to the CZ gate 
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The single-qubit measurement command corresponds to a measurement of qubit i in the basis 
|±^) = :^(|0) ± e^^|l)), with outcome si — Q associated with |+^), and outcome 1 with |— The 
measurement outcomes are usuahy referred as signals. Finahy, the corrections may be of two types, 
either Pauh X or Pauh and they may depend on any prior measurement results, denoted by 
8 — ®jeJcvSj {sj = or 1 and the summation is done modulo two). This dependency can be 
summarised as correction commands: and denoting a Pauli X and Z corrections on qubit i 
which must be applied only when the parity of the measurement outcomes on qubits j ^ J CV equals 
one (as Z^ = X^ = /). A characteristic of the MBQC model is that the choice of measurement bases 
may depend on earlier measurement outcomes. These dependent measurements can be conveniently 
written as t[M^Y^ where 

t[M!r ^ MfXfZl = Mt'^''+'\ (1) 

where it is understood that the operations are performed in the order from right to left in the 
sequence. The left (t) and right (s) dependencies of the measurement Mi are called its Z and X 
dependencies^ respectively. 

A pattern is runnable, that is, corresponds to a physically sound sequence of operations, if it 
satisfies the following requirements: (RO) no command depends on outcomes not yet measured; (Rl) 
no command acts on a qubit already measured or not yet prepared, with the obvious exception 
of the preparation commands; (R2) a qubit undergoes measurement (preparation) iff it is not an 
output (input) qubit. 

As an example, take the pattern consisting of the choices V — {1, 2}, / = {1}, O — {2} and the 
sequence of commands: 

X^^M{^Ei2Nl (2) 

This sequence of operations does the following: first it initialises the output qubit 2 in the state 
|+); then it applies /\Z on qubits 1 and 2; followed by a measurement of input qubit 1 onto the 
basis {1/V2(|0) + e-^^|l)), l/\/2(|0) - e-'^\l))}. If the result is the latter vector then the one-bit 
outcome is 5i = 1 and there is a correction on the second qubit (X2 = X2), otherwise no correction 
is necessary. A simple calculation shows that this pattern implements the unitary Jq on the state 
prepared in qubit 1, out putting the result on qubit 2, where 

^.-^(; 4:.). (3) 

The simple sequence above is a convenient building block of more complicated computations in the 
MBQC model. This is because the set of single qubit Jq (V^) together with CZ on arbitrary pairs of 
qubits can be shown to be a universal set of gates for quantum computation [2J. 

The following rewrite rules (^) put the command sequence in the standard form, where 
preparation is done first followed by the entanglement, measurements and corrections: 



E,,Xt XtZ^Ei, (4) 

tiM^yx: ^ t[Mfr+^ (6) 

i[Mf]^Z[ ^ r+t[M!r (7) 
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This procedure is called standardisation and can directly change the dependencies structure commands, 
possibly reducing the computational depth, without breaking the causality ordering ^llj . 

1.2 Determinism in MBQC 

Due to the probabilistic nature of quantum measurement, not every measurement pattern implements 
a deterministic computation - a completely positive, trace-preserving (cptp) map that sends pure 
states to pure states. We will refer to the collection of possible measurement outcomes as a branch of 
the computation. In this paper, we consider deterministic patterns which satisfies three conditions: 
(1) the probability of obtaining each branch is the same, called strong determinism] (2) for any 
measurement angle we have determinism, called uniform determinism; and (3) which are deterministic 
after each single measurement, called stepwise determinism. We will call those patterns simply 
deterministic patterns. Here since we are only working with quantum circuits we don't need to be 
concerned with other stronger notions of determinism defined for MBQC pattern such as in [2^. 

In ^] conditions over a graph (knows as flow) are presented in order to identify a dependency 
structure for the measurement sequence associated to open graph to obtain determinism. In what 
follows, we call non-input vertices as (complement of / in the graph) and non-output vertices as 

(complement of O in the graph). 

Definition 2 (Flow ll5j). We say that an open graph (G^I^O) has flow iff there exists a map 
f : and a strict partial order f over all vertices in the graph such that for all i G 

• (Fl) i <f fit); 

• (F2) if j G N{f{i)), then j — i or i j, where N{v) is the neighbourhood of v; 

• (F3) i e NUii)); 

Efficient algorithms for finding flow (if it exist) can be found in |31[ [T4]. The flow function / is a 
one-to-one function. The proof is trivial, but as this property is extensively used in this work we will 
present the proof in this paper. 

Lemma 1. Let (/, ~<f) be a flow on an open graph (G, /, O). The function f is an injective function, 
i.e. for every i G , f(i) is unique. 

Proof. Let us assume that for some i G O^, f{i) is not unique, i.e. there exists j G such that 
i ^ j but f(i) = f(j)' Then according to the flow definition: 

j G A^(/(i)) ^ i j, 
ieN{f{j))^j^fi, 

and we arrive to a contradiction because i j and j i cannot be true at the same time. Hence 
f{i) has to be unique. □ 

In the case where |/| = |0|, the flow function induces a circuit-like structure in a graph, in the 
sense that for each input qubit i G /, there exists a number n such that /^(i) G O and the vertex 
sequence 

0, /«,/[/«], /[/[/«]]},-,/"«}, (8) 

can be translated to a single wire in the circuit model. A simple example can be seen in Figure |1.2[ 
This circuit-like structure is an interesting feature of the flow function, since it allows a very simple 
translation procedure called star decomposition introduced in [15j . 
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(a) 



(b) 



Figure 1: Example of graphs satisfying the conditions of (a) flow and (b) gflow. Each vertex represent 
a qubit, where the black vertices represent measured qubits and the white ones are unmeasured 
qubits. The input qubits are represented by boxed vertices and the arrows highlight the dependency 
structure, since it points to vertices belonging to the correcting set of the measured qubit from where 
it came from. The dashed line with an arrow represent the case where the vertex from the correcting 
set is a non-neighbouring vertex (vertices not linked by an edge). 



Flow provides only a sufficient condition for determinism but one can generalise the above 
definition to obtain a condition that is both necessary and sufficient. This generalisation allows 
correcting sets with more than one element. In those cases, we say that the graph has generalised 
flow (or simply gflow). In what follows we define Odd{K) = {k , \NG{k) D K\ = 1 mod 2} to be 
the set of vertices where each element is connected with the set K by an odd number of edges. 

Definition 3 (Generalised flow [5j). We say {G^I^O) has generalised flow if there exists a map 
g : P^^ (the set of all subsets of non-input qubits) and a partial order over all vertices in 

the graph such that for all i G , 

• (Gl) if j G g{i) then i j; 

• (G2) if j G Odd(g(i)) then j — i or i ^g j; 

• (G3) i^Odd{g{i)); 

The set g{i) is often referred to as the correcting set for qubit i. It is important to note that flow 
is a special case of gflow, where g{i) contains only one element. This is a key difference regarding 
the translation of measurement patterns to quantum circuits. An example of a graph with gflow 



(but no flow) is shown in Figure 1.2 



The gflow partial oder leads to an arrangement of the vertices into layers (see below), in which all 
the corresponding measurements can be performed simultaneously. The number of layers corresponds 
to the number of parallel steps in which a computation could be flnished, known as the depth of the 
pattern. 

Definition 4 (Depth of a gfiow |14j). For a given open graph (G^I^O) and a gflow {g^-<g) of 
{GJ,0), let 



fmax^^(y(G)) ifk^O 
\max^^(y(G) \ H^kV^n) ^fk > 
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where max(X)^^ = {u ^ X s.t. Mv G X, ^{u v)} is the set of maximal elements of X according 
to ^g. The depth d^^ of the gflow is the smallest d such that V^-^ = 0; {Vk)k=o d^9 '^^ ^ partition 
ofV(G) into d^9 + i layers. 

We define the layering function of a gflow based on the above distribution of vertices into layers. 

Definition 5 (Layering function). Given a gflow {g,^g) on an open graph {G,I,0) we define 
its layering function Lg : V{G) N to be the natural number k such that i G . 

There is another useful way to understand the depth of a gflow. A gflow can be represented as a 
directed graph on top of an open graph as shown in Figure L2 The longest path from inputs to 
outputs over those directed edges corresponds to the depth of the gflow. In [14j it was shown, that a 
special type of gflow, called a maximally delayed gflow^ has minimal depth. 

Definition 6 (Maximally delayed gflow ^^). For a given open graph {G,I,0) and two given 
gflows {g, ^g) and {g\ ^^/) of (G, /, O), {g, ^g) is more delayed than {g\ ^^/) if^k, \ U^=o.../c I ^ 
I U^=o.../c ^i^^'l there exists a k such that the inequality is strict. A gflow (g, ^g) is maximally 
delayed if there exists no gflow of the same graph that is more delayed. 

We will simply refer to the maximally delayed gflow as the optimal gflow. Note that in [14j it 
was proven that the layering of the vertices imposed by an optimal gflow (^, ^g) is always unique, 
however the gflow itself might not be unique. This is an important property together with the 
following lemmas that we will exploit later for our main result on linking gflow to other known 
structures for MBQC. 

Lemma 2 (Lemma 1 from [14j). // (^, ^) is a maximally delayed gflow of (G, /, O) then = O. 

Lemma 3 (Lemma 2 from [14j). // (^, ^) is a maximally delayed gflow of (G, /, O) then (^, ^g) is a 
maximally delayed gflow of (G, /, O U V{^) where g is the restriction of g to V{G) \ (V{^ U T^q^) and 



2 Depth optimisation tools for MBQC 

Th parallel power of MBQC is proven to be equivalent to quantum circuit augmented with free 
unbounded fanout [30]. This motivates to use MBQC as an automated tool for circuit parallelisation 
as it was first presented in [TTJ. Another way to obtain parallel MBQC structure is to extract the 
entanglement graph of the pattern and obtain the optimal gflow of the graph [14j . Then one performs 
the required corrections according to this structure. Our first main result is to show the equivalence 
between these two seemingly very different technique for the patterns obtained from a quantum 
circuit, that is those with flow. More precisely we show how the effect of performing signal shifting 
optimisation (that is the core idea in [llj ) result in a maximally delayed gflow. This structural link 
shed further lights on the complicated structure of maximally delayed gflow and permit us to find a 
new efficient algorithm for finding it for the large class of patterns obtained form a circuit. 

2.1 Signal shifted flow 

We proceed with reviewing the rules for signal shifting defined in [1]: 

t[Mtr^sj[Mtr (9) 

x| sj ^ sj (10) 

Z| Sj ^ Sj (11) 
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where Sl is the signal shifting command (adding t to si) and sft/s^] denotes the substitution of si 
with t m s. Signal shifting, can be utilised to parallelise MBQC patterns and quantum circuits [iT]. 
The rest of this section is focused on various structural properties of the signal shifting. 

As can be seen from the above rules, signal shifting rewrites the X- and Z-corrections of a 
measurement pattern in a well defined manner. In particular, it will move all the Z-corrections 
to the end of the pattern, thereby introducing new X-corrections when Rule [TO] is applied. It is 
proven in [Tlj that signal shifting will never increase the depth of an MBQC pattern, although it 
can decrease it. In the case when the depth decreases, it is the consequence of the removal of the 
Z-corrections on the measured qubits by applying Rule |9| 

The Rules [9] - [TT] can be interpreted in the following way. Signal shifting takes a signal from 
a Z-correction on a measured qubit i (Rule ^ and adds it to the corrections that depend on the 



outcome of the measurement of i (Rules [1Q| - 11). When the signal moves to an X-correction 



command, then it won't propagate any further. If the signal was added to another Z-correction 
of a measured vertex, then signal shifting can be applied again until no Z-corrections are left 
on non-output vertices. Therefore signals move along a path created by the Z-corrections. The 
propagation of signals in an MBQC pattern can be described by a Z-path as defined below. 

Definition 7. Let M be a measurement pattern on an open graph {G,I,0). Then we define a 
directed acyclic graph, called Gz, on the vertices of G such that there exists a directed edge from i to 
j iff there exists a correction command Z|' in M. A path in Gz between two vertices v and u is 
called a Z-path. 

The above definition allows us to state a simple observation about connectivity of a graph with 
flow. 

Lemma 4. // (/, ~<f) is a flow on an open graph (G, /, O), and there exists a Z-path from vertex i 
to vertex j, then the vertices i and f{j) cannot be connected. 

Proof. The existence of a Z-path from i to j implies that i j. The Z dependency graph is an 
acyclic graph, thus i ^ j. If i would be connected to /(j), then according to the flow property (F2): 

i^N{f{j)) A i^j ^ j^fi 

Now we have two contradicting strict partial order relations i j and j i. Therefore i cannot 
be connected to /(j). □ 

Recall that the addition of signals is done modulo 2, therefore, if an even number of signals from 
a measured vertex i is added to a correction command on vertex j, the signals will cancel out (since 
Z^ = = /). Furthermore, it is evident from the rewrite Rules of [o]- [ll] that after signal shifting, 
the measurement result of vertex i will create a new X-correction over vertex j if there exist an odd 
number of Z-paths from z to a vertex k that j is X-dependent on directly before signal shifting. 
Similarly a new Z-correction from i to j will be created if there exists an odd number of Z-paths 
from i to j. Either way, the number of Z-paths from a vertex i to another vertex j, denoted as QU)^ 
can be used to determine if the signal from i should be added to a correction. We define (i{i) to be 
1 to simplify further calculations and definitions in this paper. The importance of the number of 
Z-paths will manifest itself in the next subsection, when the relation between signal shifting and 
gflows is studied. 

We define a new structure called the signal shifted flow (SSF), and show that it satisfies the 
three gflow properties in Definition [3| Before constructing the SSF, some definitions and lemmas are 
needed to justify our definition. Note that if an open graph (G, /, O) has a flow (/, ^/), then we 
can write the MBQC pattern of a deterministic computations on this open graph as [15j : 
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n {^h^Nifmm^r) ^oNjc (12) 



12 



we 



where the product follows the strict partial order ^/ of the flow (/, -<f). From Equation 
see that a Z-correction on a vertex j depending on the measurement outcome of another vertex i 
appears only if j is a neighbour of f{i). This is formally stated in the next corollary as we will refer 
to it several times. 

Corollary 1. // (G, /, O) is an open graph with a flow (/, ^f), then there exists a Z- correction from 
vertex i to another vertex j iff j G N(f(i)) \ {i}. 

We define Z -dependency neighbourhood of a vertex j to be the set of vertices from which j is 
receiving a Z-correction from. This set has an explicit form given as Nz{j) = {A: E 0^\f{k) G 
N {j)\{f (j)}} ^ this is due to the following facts: for all vertices k G O^, from fiow definition f{k) exists 
also since f{k) ^ {/(j)} hence k cannot be equal to j and moreover since f{k) G N{j) ^ j E N{f{k)) 
therefore according to Corollary [T] there exists a Z-correction from k to j. It is easy to see, that 
Ci{j) can be written as: 

keNzU) 

There exists a Z-correction from every k G Nz{j) to j. These Z-corrections can be used to extend 
every such Z-path to k to reach j. If i is in the sum, then because (i{i) = 1 the correct number of 



Z-paths is obtained with Equation 13 



We now present the complete algorithm (Algorithm [T]) for signal shifting a fiow pattern shown in 
Equation [T2j We keep in mind that the order in which we apply the signal shifting rules does not 
matter [15j. 



as 



Proposition 1. Given the measurement pattern P of a flow (/, as defined in Equation 12 
input to Algorithm^ the output will be the signal shifted measurement pattern of P. 

Proof We will prove this proposition by showing that: 

• Algorithm [T] terminates. 

• Every step in Algorithm [T] that modifies the pattern P^ is a valid application of a signal shifting 
rewrite rule. 

• The output of Algorithm [H the pattern P\ is signal shifted. 

We begin by showing that Algorithm [T] will terminate. The first "while" loop will obviously 
terminate, as we decrease the number of elements on each loop iteration and never add anything 
to the set B. The second "while" loop will not terminate only if some Z^"- command will be added 
to the pattern an infinite number of times. As the underlying graph is finite and a Z^^ command 
represents a directed edge in the Z-correction graph, this implies the existence of a cycle in the 
graph, however this is impossible according to the fiow definition. The "for" loop in the algorithm 
terminates because the graph itself is finite, hence Algorithm [l] has to terminate. 

For Algorithm [T] to actually perform the signal shifting, its operations have to be either trivial 
commuting rules or the three signal shifting Rules [9] - [TOj As can be easily seen from the algorithm, 
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Algorithm 1: SignalShift 



Input: A measurement pattern P with flow (/, ^f) as deflned in Equation 

Output: The signal shifted pattern of P. 

begin 

P' = P; 

while B ^ do 

select any vertex i G B which is smahest according to -</ ; 

B = B\{i}; 

while 3k e B s.t Zf' G do 
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Move the Zf 



command next to the M^^ command; 



Use Rule [o] on P' to create the signal command Sj^ , 
// Removes the Z^'' command from P'; 
Use Rule 10 on P' to create a new A^'^^^J command; 



for j G N{f{k)} \ {k} do 
^ Use Rule 
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on P' to create a new 



command. 



Move Si' 



to the end of the pattern and remove it. 



the operations done are indeed the signal shifting rewrite Rules |9] - [T0| We still need to prove, that 
these rules can be applied in the order shown in the algorithm. Obviously we can use Rule |9] on line 
8 to create the signal command due to the fact that A: E P C and that every non-output qubit is 
measured. Hence we have the measurement required for the creation of the signal command in the 
pattern. We know that has to be in the pattern after the command M^^' and before M^^ . The 
entanglement and creation commands are the first commands in the pattern and we do not need to 
move the command past them. Hence we only need to move past measurement commands 
on qubits that are not i and k and other correction commands. These can be done trivially and 
hence we can always move the command next to M^^ to apply Rule [oj 

Next we want to move the newly created SI command to the end of the measurement pattern. 
To do that we need to commute it past the commands that appear after it. The only commands S^ 
commutes non-trivially with are the ones that depend on the measurement of qubit k as can be seen 
from Rules [9] - [TT] Those are the X- and Z-corrections depending on the measurement outcome of 



qubit k. According to Equation 12 there is exactly one such A-correction in the pattern P, namely 
A^^^^J . Also the previous steps of the algorithm could not have created any dependencies from qubit 
k. The Z-correction commands have only been created depending on vertices that we already moved 
from B. Therefore we need to create exactly one new A-correction command using RulelTOl We also 



look at the Z-corrections depending on k and from Equation 12 we see that in the original pattern 
these are on vertices from the set N{f{k)) \ {k}. As for the A-corrections we also have not created 
any new Z-corrections from k in the previous steps of the algorithm. Hence this is exactly the set of 



corrections we need to commute with and apply Rule 11 We are only left with commands after SI 
in the pattern that commute trivially with S^. We can move the command at the end of the pattern. 
The signal command at the end of the pattern does not influence the computation and we will not 
add any new commands to the end of the pattern. Hence we can remove the S^ command from the 
pattern. 

Finally we show that no more signal shifting rules can be applied after the completion of Algorithm 
[H i.e. the pattern P' is signal shifted. We eliminate all Z-corrections acting on a non-output qubit 
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depending on a vertex i after removing it from the set B and will afterwards never create any new 
Z-corrections depending on that vertex. At the end of the algorithm the set B is empty, hence there 
cannot exist any non-output qubit that has a Z-correction acting on it and Rule [9] cannot be applied 
anymore. Moreover, since every signal command is at the end of the pattern, we cannot apply the 



Rules 10 and 11 neither, that completes the proof. □ 



We consider any trivial commutation of a pattern commands resulting to an equivalent pattern. 
Therefore the above algorithm defines the unique pattern obtained after signal shifting. Note that 
Algorithm [T] works almost like a directed graph traversal, where there is a directed edge from vertex 
i to A; iff there exists the command in the measurement pattern. The only difference from a 
classical directed graph traversal is that we allow visiting of a vertex more than only once. Hence we 
will traverse through every different path in the graph however we do that exactly once. 

As mentioned before, the evenness of the number of Z-paths can be used to determine if a signal 
is added to a correction command. Let parity{n) be the function that determines the oddness or 
evenness of the integer n, i.e. parity (n) = n mod 2. Then if an open graph has a flow, the oddness 
of CiU) t)e found as described in the following lemma. 

Lemma 5. For every two vertices i and j in an open graph (G^I^O) with flow (/, ^f) 

parity{Q{j)) = \k e {Nz(j)\parity(Q(k)) = 1}| mod 2 

i.e. parity{Q{j)) depends only on the number of vertices in the Z-dependency neighbourhood which 
have odd number of Z -paths from i. 

Proof. The oddness of CiU) can be written as 

parity{Ci{j)) = ^ Ci{k) mod 2 = 
\keNz{j) J 

^ {Ci{k) mod 2) mod 2 
keNzU) 

^ {Ci{k) mod 2) mod 2 

{keNz(j)\l=Q(k) mod 2} 

= \{keNz{j)\l = Q{k) mod2}| mod 2 = 
-\{keNz{j)\parity{Q{k))^l}\ mod 2 

□ 

All these notions will allow us to define the structure of the pattern after signal shifting is being 
performed. 

Proposition 2. Given a flow (/, ^ /) on an open graph (G, /, O), let s be a function from ^ P^^ 
such that j G s{i) iff parity {(i{f~^{j))) = 1. Also define Lg to be a layering function from V{G) 
into a natural number: 

L,(z) = VzgO 
Ls{i) = max(L5(j) + 1) Vz ^ O 



12 



Define the strict partial order -<s with: 



i^sj ^ Ls{i)>Ls{j) 

Then, the application of signal shifting Rules^ -[77] over an MB QC pattern with flow (/, ^j) will 
lead to the following pattern: 

p = n n (x'(V)Mf') EgNjc (14) 

Proof. The proof is divided into three parts. First we wih show that signal shifting creates exactly 



the pattern commands shown in Equation 14 We proceed by showing, that the layering function Lg 
is defined for every i G V{G). Lastly, we need to prove that using the partial order derived from 
Ls for ordering the commands as in Equation [M] gives a valid measurement pattern. 

Note that the preparation commands {Nf)^ entanglement commands {Eg) and measurement 



commands {M^"-) are the same for Equations 12 and 14 Because signal shifting would not change 



these commands (Rules [9|- [Tlj) these are as required for a signal shifted pattern. Hence we need 
only to consider the correction commands. 

We will look at the correction commands that would appear in a signal shifted pattern. We do 
this by examining the signal shifting algorithm (Algorithm [T]) . As mentioned before, the algorithm 
works as a directed graph traversal, in a way that every distinct path is traversed o. As seen in the 
algorithm every Z^* correction acting on a non-output qubit is removed from the pattern. This is in 
accordance with the proposed pattern in Equation [T4} Let us examine which new corrections are 
created. 

The number of newly created depends on the number of times we enter the first loop with 



command Z^j^_-^^.y As the algorithm is a directed graph traversal algorithm, this happens as many 

times as there are different paths over the Z-dependency graph from i to f~^{j). Because the same 
two Xi'' corrections cancel each other, hence a new X-correction appears in a signal shifted pattern 
only if parity {Q{f~^ (j))) = 1. We also note that no new X(/(z))^^ correction is created since there 
exist no Z-path between i and f~^{f{i))- On the other hand Algorithm [l] leaves the already existed 
X corrections unchanged and moreover since we have defined Q{i) = 1 therefore f(i) G s{i). This 
implies that the set s{i) does indeed contain all the vertices that have an X-correction depending on 
Si after signal shifting is performed. 

The number of newly created Z-corrections on an output vertex j depending on a vertex i 
appearing in the signal shifted pattern is equal to the number of different paths from i to j. The 
difference with non-output qubits is that these will not be removed through the process of signal 
shifting. As with X-corrections, two Z-correction commands on the same qubit will cancel each 
other out and hence the existence of a in the final pattern depends on the parity of the number 
of paths from i to j. This can be written in short as: 

^Si-parity{CiU)) 



Hence the measurement pattern in Equation 14 has exactly the same commands as the signal shifted 



pattern in Equation 12 



Another thing we need to proof is that the layering function Lg is defined for every i G V{G). 
As proven above, the X-corrections depending on the measurement of qubit v correspond to the set 
s{v). Hence we can interpret the definition of Ls{v) as finding the maximum value of Ls for every 
vertex that has an X-correction from v and adding 1 to it. The recursive definition of Ls{v) is well 
defined, if for every non-output qubit we can find a path over X-corrections ending at an output 
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qubit. We know that signal shifting of a vahd pattern creates another vahd pattern. This imphes 
that the X-corrections cannot create a cychc dependency structure and hence every path over the 
X-corrections has an endpoint. Moreover such a path cannot end on a non-output qubit k since 
f{k) G s{k) and one could always extend that path with f(k). Therefore Ls{v) is well defined. 

Finally, it is easy to show that the partial order as used in Equation [14] gives a valid ordering 
of the commands. Every vertex j that has an X-correction depending on the measurement of qubit 
i has a smaller Lg number and hence i j- This way no X-correction command acts on an already 
measured qubit and because the Z-corrections are applied only on output qubits, the correction 
ordering is valid. Every other command is applied before the measurement command and hence the 
pattern in Equation [M] is a valid measurement pattern. □ 

Given an open graph with a flow, we refer to the construction of the above proposition as its 
corresponding signal shifted flow (SSF). The main theorem of this section states that every SSF is 
actually a special case of a gflow. 

Corollary 2. // (G, /, O) is an open graph with flow (/, ^f) and SSF (5, ^5) then for every vertex 
i and j such that f(j) G s{i) \ {/(i)}; we can find another vertex k, such that f{k) G s{i) H N(j). 

Proof If f{j) G then from the Proposition [2] of SSF we can conclude that parity {^{j)) — 1. 
We know that j ^ i from the assumptions. Lemma [5] says, that there must exist at least one other 
vertex k from which j has a Z-correction, such that parity{Q{k)) = 1. The flow definition says that 
j must therefore be a neighbour of f{k). Definition [2] of SSF states that f{k) must therefore be in 
hence f{k) G s{i) H N{j). □ 

Theorem 1. Given any open graph (G^I^O) with flow (/, ~<f), the corresponding signal shifted flow 
(5, ^5) is a gflow. 

The proof is based on the following lemmas, demonstrating that 5 is a gflow by satisfying all the 
properties of Definition [sj The first property of gflow (property Gl) is satisfied by SFF implicitly 
from Definition [2j i.e. for every i G V{G) it holds that i j if j ^ s{i). Consider the second gflow 
property (G2), i.e. if j G Odd{g{i)) then j — i or i j- We will show that every vertex with odd 
many connections to s{i) has to be either i itself or an output qubit. This is a stronger condition 
than is needed to show G2, but as shown in Section [3} necessary for creating compact circuits from 
SSF. 

Lemma 6. // (5, ~<s) is an SSF then every non-output vertex v ^ i connected to s{i) has an even 
number of connections to s{i), i.e., 

\/v G N{s{i)) \0 A v^i ^ V ^ Odd{s{i)) 

Proof. Let 7^ z be a vertex connected to we show the following two sets have the same number 
of elements. 

{k G Nz{v) I parity{C^{k)) = 1} and s{i) H N{v) \ {f{v)} 

For every j G s{i) H N{v) \ {/('y)}, we prove f~^{j) is the unique element in 

{k^Nz{v) \par^ty{C^{k))^l} 

Because j G s{i) from Proposition [2] there must exist f~^{j). Also since j G N{v) therefore 
V G N{j) — N{f{f~^{j))). Moreover since j 7^ /('^), Corollary [l] implies the existence of a Z- 
correction from f~^{j) to i.e. f~^{j) G Nz{v). Proposition [2] says that because j G it must 
hold that parity {Ci(f-^{j))) = 1. Therefore f~^{j) G {A: G Nz{v) \ parity{Ci{k)) = 1}. 
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On the other hand, for every vertex u G {k ^ Nz{v) \ parity{(i{k)) = 1}, as parity{(i{u)) = 1 
then from Proposition [2] we haye f{u) E s{i). Also f{u) G N{v) because of Corohary [l] and finahy, 
f{u) ^ f{v) because v cannot have a Z-correction from itself, i.e. v ^ Nz- Hence it holds that 
f \u) E s{i) n N{v) \ {f{v)}. Therefore 

\s{i) n N{v) \ {f{v)}\ = G Nz{v) \ parityiQik)) = 1}| 

According to Lemma [5] par it?/ (("^(i;)) = \{k G Nz{v) \ parity{Q{k)) = 1}| mod 2. If parity {Q{v)) = 
then s{i)nN{v)\{f{v)} must have an even number of elements. Proposition [2] says that f{v) cannot 
be in s{i) and therefore v can have only even number of connections to s{i). If parity {Q{v)) = 1 
then we know that s{i) H N{v) \ {f{v)} must have an odd number of elements. If f{v) exists 
it must be in s{i) because of Proposition [2j In the case of f{v) G we can conclude that 

parity{\s{i) H N{v)\) — and v has even many connections to s{i). On the other hand if f{v) does 
not exist, v has to be an output qubit because the flow function / is defined for every non-output 
vertex. The only possibility of k having odd many connections to s{i) is therefore if k is an output 
vertex, which proves the lemma. □ 

The next lemma directly proves that an SSF also satisfies the last gflow property (G3) which 
states that i G Odd{s{i)). 

Lemma 7. // (5, ^5) is an SSF, then for every i G it holds that i G Odd{s{i)). 

Proof. First we show that, performing signal shifting creates new X-corrections only between 
unconnected vertices. Recall that signal shifting creates a new X-correction between vertices i and j 
iff there exists a Z-path from i to f~^{j) and an X correction from f~^{j) to j therefore from the 
Flow definition we have: 

Let us assume that there exists an edge between i and j. According to the Flow definition we have 
that 



This contradicts the partial order i f~^{j) j of the Flow (/, ^f) and therefore there cannot 
be an edge between vertices i and j. 

Next we claim that there is exactly one edge between i and s{i). According to Definition [2] of 
SSF, the set s{i) consists only of the vertex f{i) and the vertices to which signal shifting created a 
new X dependency from i. We showed that signal shifting does not create X dependencies between 
connected edges. Hence, f(i) is the only vertex in s(i) that can be connected to z, and there must 
be an edge between i and f{i) because of the flow property (F3) {i G N{f{i))). □ 

Proof. Now to obtain the proof of Theorem [TJ we note that the definition of SSF implies the gflow 
property (1.). Lemma |6] implies that every SSF satisfies the gflow condition (2.). As the third and 
last gflow condition is satisfied by SSF according to Lemma [7| SSF is indeed a gflow and Theorem [T] 
holds. □ 

The above theorem for the first time presents an structural link between two seemingly different 
approach for parallelisation, gflow and signal shifting, for those patterns having already flow. As 
mentioned in the introduction this is the key step in obtaining our simultaneous depth and space 
optimisation. The next section explores further the link with gflow, showing optimality of SSF in 
parallelisation. 



15 



2.2 Influencing paths 

The notions of influencing walks and partial influencing walks on open graphs with flow was 
introduced in [11] to describe the set of all vertices that a measurement depends on. An influencing 
walk starts with an input and ends with an output vertex, a partial influencing walk starts with an 
input vertex but can end with a non-output vertex. We will use a modified definition of influencing 
walks that can start from any non-output vertex i and end at any vertex j G s{i) and call it a 
stepwise influencing path. This will allow us to conveniently explore the dependency structure of a 
pattern with SSF. 

Definition 8. Let (5, ^5) be an SSF that is obtained from a flow (/, -<f) of an open graph (G, /, O) 
and vertices i and j in V{G) such that j G s{i). We say that a path between vertices i and j is an 
stepwise influencing path^ noted as pi{j), iff 

• The path is over the edges of G. 

• The first two elements on the path are i and f{i). 

• Every even-placed vertex k on the path pi{j), starting from f{i), is in s{i). 

• Every odd-placed vertex on the path pi{j) is the unique vertex f~^{k) of some k G s{i) such 
that k is the next vertex on the path pi{j)- 

It is easy to see that every second edge, in particular the edges between f~^{k) and k G s(z), in 
the stepwise influencing path is a flow edge. Hence the path contains no consecutive non-flow edges. 
If we restrict the first vertices of the stepwise influencing path to be input vertices, the stepwise 
influencing path would be a partial influencing path, but not vice versa. Stepwise influencing paths 
are useful because of their appearance in the SSF as proven by the following lemma. 

Lemma 8. Let (5,^5) be an SSF obtained from a flow {f,^f) of an open graph {G^L,0) and 
vertices i and j in V{G) such that j G s{i). Then there always exists a stepwise influencing path 

piU)- 

Proof We start by constructing such a path backward from j to i. We select j and f~^{j) as the 
last two vertices on the path and apply Corollary |2] to find the vertices on the path, until we reach i. 
The formation of cycles is impossible, as this would imply a cyclic dependency structure, impossible 
for a flow. We have to reach i as the set of vertices we choose from is finite. □ 

The above lemma will be used in Section [3] to obtain compact circuits from SSF. Note that 
there might be more than one stepwise influencing path from i to j. We conclude the section about 
influencing paths with the following two lemmas which will be used to prove the optimality of SSF. 
First, the structure of stepwise influencing paths imposes a strict restriction on the way a vertex on 
the stepwise influencing path can be connected. 

Lemma 9. Let pi{j) be a stepwise influencing path from i to j in an open graph (G, /, O) with flow 
(/, ^f) and corresponding SSF (5, ^5). Then f~^{j) is the only odd-placed vertex in pi{j) that j is 
connected to. 

Proof. According to the definition of stepwise influencing path, for every three consecutive vertices vi^ 
V2-, '^a in pi{j) such that vi and 173 are odd-placed we have that V2 — f{vi) and vs G N{v2) = N{f{vi)). 
According to Corollary [T] there must exist a Z-correction from vi to vs. Therefore the odd-placed 
vertices in pi{j) are on a Z-path from i to f~^{j) and obviously from every odd-placed vertex in 
pi{j) there exists a Z-path to f~^{j). Lemma [4] says that j cannot be connected to any of the 
odd-placed vertices in pi(j). □ 
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The previous lemma shows, that the stepwise influencing paths can be used to describe some 
properties of the connectivity in open graphs with SSF. The next lemma (illustrated in Figure [2]) 
will explain how a stepwise influencing path can be extended. 

Lemma 10. Let {G^I^O) he an open graph with flow (/, ^j) and corresponding SSF (5, ^5) and let 
i and j be two non-output vertices of the open graph such that f{j) G s{i). If v ^ ^ s{i) \ {f{j)} 
then every stepwise influencing path pi{v) can be extended by the vertices j and f{j) to create another 
stepwise influencing path pi{f{j))- 

Proof Adding j and f{j) to pi{v) satisfies the conditions for stepwise infiuencing paths. There 
exists an edge between vertices j and v and vertices j and /(j), hence it is a valid path. Moreover, 
f{j) G s{i) would be an even-placed vertex on the extended path, and j would be the unique 
oddly-placed vertex with f{j) E s{i). □ 




Figure 2: Extending a stepwise infiuencing 
path ending at vertex v according to Lemma 

m 



Figure 3: For every strict subset S of 
s{i) containing f{i) we can find a vertex v 
in the odd neighbourhood of S such that 
f{v) is not contained in S. This is proven 
in Lemma [TTl 



2.3 Optimality of signal shifting 

Given an MBQC pattern with gfiow, finding the maximally delayed gfiow of its underlying graph 
could potentially further reduce the depth of the computation [14j. A natural question that arises is 
how SSF is linked with the optimal gfiow. In this section, we prove that if the input and output sizes 
of the pattern are equal, then SSF is indeed the optimal gfiow. Hence we can conclude the most 
optimal parallelisation that one could obtain via translation of a quantum circuit into an MBQC 
pattern is achieved by the simple rewriting rules of SSF. This will also lead to a more efiicient 
algorithm than the one presented in [14j for finding the maximally delayed gfiow of a graph as we 
discuss later. 

Theorem 2. Let {G,I,0) be an open graph with flow (/, such that \I\ = |0|. Let (5, ^5) be the 
SSF obtained from (/, ^/). Then (5, ^5) is the optimal gflow of (G, /, O). 
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The proof of the theorem is rather long, an outhne is presented below. A general reader could 
omit the next subsections, however various novel constructions has been introduced in the proof 
that could be explored for other MBQC results and hence could be valuable for an MBQC expert. 
In Section 2.3. 1| we show that the penultimate layers of an optimal gflow and an SSF of an open 
graph where |/| = |0|, are equal. Next we introduce the concept of a reduced open graph in Section 



2.3.2[ We prove two key properties of the optimal gflow and SSF of the reduced open graph. This 



highlights the recursive structures of the gfl ow and SSF leading to the possibility of extending these 
notions to new domains In Section 2.3.3 we put the pieces together, by showing that the previous 
properties imply that reduced gflow (implicitly also optimal gflow and SSF) layers are equal to the 
original gflow layers from layer 1 onward. This allows us to construct a recursive proof for Theorem 
|2j which we present in Section 2.3.4 



2.3.1 The last two layers 

The equality of the last layers of an SSF and optimal gflow follows from Lemma |2] and Proposition [2] 
- the last layer of an optimal gflow and an SSF is always the set of output vertices. What is left to 
prove is that the penultimate layers are also equal, for doing so we need the following properties 
of open graphs with SSF. An illustration of the property proven in the flrst of the two lemmas is 
shown in Figure |3| 

Lemma 11. Let {G,I,0) be an open graph with flow {f,^f) and corresponding SSF (5,^5). // 
i G then for every strict subset S of s{i) containing f{i) there must exist a non-output vertex v 
that is oddly connected to S such that f{v) G s{i) \ S, i.e. 

Vz G V5 C s.t. f{i) eS 3ve Odd{S) s.t. f{v) G s{i) \ S 

Proof. If s{i) — {/(i)} the lemma holds trivially, as there does not exist any nonempty strict subsets 
of s{i). Consider the case where s{i) contains more than one element and 5 is a strict subset of s{i). 
Then we select any vertex j ^ S from s{i) and look at the stepwise influencing paths from i to j. 
Note that there might be more than one such path. We move backwards from j towards i over the 
stepwise influencing paths in the following way: 

1. Move by two vertices 

1.1 If possible, choose any stepwise influencing path where the previous even-placed element 
is not in S and move to that element. 

1.2 If the previous even-placed elements in all the stepwise influencing paths from i to j 
are in 5, then stop. 

2. Repeat step 1. 

Let u be the vertex to where we moved using the above process, u has to exist because of the way 
we initially selected j. There are a couple of other observations that we can make about u. First, 
u G s{i) \ 5, because of the selection of j and the way we moved on the paths. Second, u cannot be 
the flrst even placed vertex on a stepwise influencing path from i to u because the flrst element is 
f{i) G S (according to Deflnition [s]). Third, for every stepwise influencing path ending in the 
previous even-placed vertex has to be in S as otherwise we could have moved one more step towards 
i. 

^For example, the authors are currently exploring this structure to define the concept of partial flow, for patterns 
with no deterministic computation. 
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Considering the previous three observations we can show that the vertex v — f~^{u) must be oddly 
connected to S. We begin by noting that v cannot be connected to any vertex k G s{i) \{SU {f{v)}. 
Otherwise, according to Lemma [TOl we could extend any stepwise influencing path ending at k with 
V and f{v). Hence k ^ S U {f{v)} would then be an even-placed vertex on a stepwise influencing 
path from i to f{v). In particular, k would be the second to last even-placed vertex on a stepwise 
influencing path from i to f{v) = u Every such vertex, except f{v) itself, is in S as mentioned before. 
Because, according to Lemma [oj v has to be evenly connected to it has to be oddly connected 
to S and Lemma [TT] holds. □ 

Next we need to show that every non-input vertex i has a corresponding unique vertex f~^(i)^ 
this is only true for those graphs with |/| = |0|. 

Lemma 12. // (/, ^f) is a flow on an open graph (GJ,0), then \I\ — \0\ iff for every j G there 
exists f~^{j). 

Proof. First, if |/| = \0\ then also — \0^\. The flow deflnition uniquely deflnes f{i) for every 
i G and therefore f~^{j) is uniquely deflned for some, but not necessarily for all, vertices j G . 
The number of vertices for which / is deflned must equal the number of vertices for which f~^ is 
deflned and because = f~^ must be deflned for every element in 

Second, Let us consider the case when for every j G there exists f~^{j). The number of 
elements for which is deflned equals the number of elements / is deflned for. / is by Deflnition 
deflned for every element in . Hence — \0^\ which implies that |/| = |0|. □ 



Note that the above requirement, i.e. the existence of f~^{i)^ is the only reason why our proof of 
Theorem [2] fails if |/| 7^ \0\. We conjecture that by padding the input with necessary ancilla qubits 
without changing the underlying computation we could extend the above theorem to the general 
graphs. However the proof of such result is outside of the scope of this paper and not relevant for 
the optimisation of quantum circuit. 

Note that because of Deflnition [6] if a gflow is not optimal, its penultimate layer has to either be 
equal to the penultimate layer of the optimal gflow or there exists a vertex in the penultimate layer 
of optimal gflow that is not included in the penultimate layer of the other gflow. In the proof of the 
main result we assume that the penultimate layers are not equal, hence we could choose a vertex 
with particular properties (described in the next two lemmas) to derive a contradiction. 
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Lemma 13. Let (G^I^O) be an open graph where \I\ = \0\ with flow {f^^f), corresponding SSF 



(5, ^5) and a gflow 



such that Vq ^ 



O. Assume there exists a vertex i ^V^ ^ \ then 
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• 9{i) C O 

• gii)ns{i) c s{i) 

• fii) e gii) 

and there exists a vertex bo such that 

• bo e Odd{g{i) n s{i)) 

• f{bo)(Esi{)\gi{) 

Proof. Because i is in V^^ the set g{i) must be a subset of V^^ = O according to Definition ^ 
Proposition [2] implies that V^^^ = O. This and the fact that i ^ V{^^ imphes that s{i) is not a subset 
of the output vertices O = V'o^^. Therefore there must exist a non-output vertex in s{i) and, because 
g(i) C O, this vertex cannot be contained in g(i). Thus the intersection of s{i) and g{i) cannot be 
equal to s(i) and g{i) H s(i) C s(i). 

We now show that f{i) G ^(i). Let us assume that f{i) ^ g{i), and choose a vertex ai G ^(i) 
connected to z, such a vertex has to exist because the gflow definition says that i is oddly connected 
to g{i). As ai G g{i) then by the gflow definition ai cannot be an input qubit. According to Lemma 
12, there must exist a vertex f~^{ai) to which ai is connected to. By the definition of flow, f~^{ai) 



cannot be an output vertex and thus is not in layer V^^ . As g{i) C O this also means f~^{ai) ^ g{i). 
On the other hand f~^{ai) is connected to ai G g{i). Because i G V^^ and f~^{ai) ^ V^^ we know 
from Definition [4] that i y^g f~^{ai). As f~^{ai) is connected to g{i) we can conclude from the 
gflow definition that f~^{ai) has to be evenly connected to g{i) and therefore has at least one more 
connection to a vertex a2 G g{i). 

Using the same argument for a2 as for ai we can say that there must exist f~^{a2) ^ g{i) to 
which a2 is connected to. Let us assume that f~^{a2) is not connected to ai. This means it has 
only one connection to the set A2 = {ai, a2} ^ g[i) and is therefore oddly connected to it. We can 
continue this procedure of selecting vertices from g{i) until we select a vertex such that f~^{an) 
is connected to at least one vertex aj in = {ai, . . . a^_i} C g(i). If this happens we can no 

longer say with certainty that f~^{an) is oddly connected to ^ ^(z), which means we cannot 
select any more elements from g(i) using this method. Because (G, /, O) is a finite open graph we 
must find this an in finite number of steps. 
We created the set An in such a way that: 

Vj G {1, 2, . . . , n - 1} f-\aj) e N{aj+i) = Nifif-\aj+i))) 

Hence we have a Z-correction from every f~^{aj-^i) to f~^{aj) and thus there exists a Z-path 
from f~^{an) to every f~^{aj) such that aj G An-i and , because of Lemma [sj f~^{an) cannot 
be connected to any vertex in An-i- This leads to a contradiction with the assumption that it is 
connected to at least one vertex in An-i. Therefore our initial assumption that f{i) ^ g{i) must be 
false and g{i) must contain f(i). 

From the definition of SSF we have that f{i) G s{i) and therefore also f{i) G g{i) H s{i). Now 
we know that g{i) H s{i) is a strict subset of s{i) containing /(i); the existence of 60 follows from 
Lemma [TTl □ 



Now we prove that if we have a vertex with the same properties as bo in Lemma 13 and a 
(possibly empty) subset A of vertices with particular properties (which will be defined in the next 
lemma) we can always increase the size of A and find another vertex with properties of 60 • This 
would imply the possibility of increasing the size of A to infinity and will give us the contradiction 
we need. 



20 



Lemma 14. Let {G,I,0) be an open graph where \I\ — \0\ with flow [f^^f), corresponding SSF 
(5, ~<s) and a gflow ^g). If we have a vertex i in the open graph such that 

9ii) C O 

g{t)r\s{i) c s{i) 
m e g{i) 

and if we have a subset A C g(^i) and another vertex bo such that 
bo G Odd{g{i) D s{i)) 
f{bo)es^)\g{i) 



then 



there exists another vertex Cq and a non empty set B C g{i) such that 
BnA = $ 

Co G Odd{g{i) n s{i)) 

f{co) e s{i) \ g^) 

yjeAuB 3 c^^f-\j) 

Proof The proof consists of three steps: we start by constructing the set B] we proceed with finding 
the vertex cq; and finally we prove that cq has the required properties. 

Define S = g{i)ns{i)^ since f{bo) exists hence bo cannot be an output vertex. Also since g{i) C O 
therefore bo is not in g{i). As bo ^ O — V^^ and g{i) C O we can conclude from Definition [4] 
that i G V^^ and i -/<g bo. Therefore according to the gflow definition, bo must be in the even 
neighbourhood of g{i). We also know from the initial conditions of this lemma that bo is in the odd 
neighbourhood of g{i) H s{i). Thus there has to exist a vertex vi in g{i) to which bo is connected to, 
but which is not included in g{i) D i.e. vi G g{i) \ s{i). As g : O P , i^i G g{i) cannot 



be an input qubit and because /~ exists for every non-input qubit according to Lemma [12j there 
must exist a vertex f~^{vi) — bi. It is also important for the later part of the proof to note that 
f{bi) — vi ^ A. This is due to Lemma |4j which implies that 60 cannot be connected to any vertex 
in A. 

Define Bo — S and consider the case when bi is evenly connected to Bo- Remember that the flow 
property (F3) says that there is always an edge between bi and f{bi). This means that bi is oddly 
connected to Bi — U Bo which is a subset of g{i). But again because of the gflow property 

(G2) we have that bi must be evenly connected to g{i). Thus there must exist another vertex 62 such 
that 61 is connected to f(b2) G g(i) otherwise 61 could not be in the even neighbourhood of g{i). 
If 62 is evenly connected to Si, it must be oddly connected to B2 — {/(&2)} U Bi which is again a 
subset of g{i). If 62 is oddly connected to B2 there must exist a vertex 63 such that 63 is connected 
to f{b'i) G g{i) \ S2, otherwise 62 could not be in the even neighbourhood of g{i). We can continue 
this scheme until we get to vertex bn that is oddly connected to Bn-i- As Bn = {/(^n)} U ^n-i and 
there exists an edge between bn and f{bn) we get that bn must be evenly connected to Bn- Such 
vertex bn must exist, otherwise we could continue selecting elements from g{i) infinitely, but (G, /, O) 
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is a finite open graph. We select B = Bn \ S. Recall that f{bi) must exist, therefore B must have at 
least on element. 

Next we show bn is oddly connected to S. We note that we have the following: 

Vie{l,2,...,n} bjeNifibj)) A bj.ieNifibj)) 

Corollary [l] implies that for every j > there exists a Z-correction from bj to bj-i. Thus we have 
a Z-path from b^ to every other bj where j < n, hence from Lemma [4] we conclude bn cannot be 
connected to any vertex f{bj) G Bn-i where j < n. The number of edges that connect the vertices in 
Bn-i to vertex bn has to be the same as the number of edges between vertices of S and 6^, because 
Bn-i — {/(^i), /(^2), • • • , f{bn-i)} U S. As bn was oddly connected to S^-i, it must also be oddly 
connected to S. Note that however bn does not have the required properties for cq, but will be used 
to find such a vertex. 

The gflow definition says that bn must be evenly connected to s{i). It is also oddly connected to 
s{i) n g(i) hence there must exist a vertex c G s(i) \ g(i) to which bn is connected to. According to 
Lemma [s] there exists a stepwise influencing path pi{c) and due to Definition [sj f{i) has to be on on 
this path. Therefore there exists at least one element in pi{c) that is in S. Let /(ag) be the last 
element of the path pi{c) in S. 

Define ai to be the vertex in pi{c) that comes after /(ao). We know that ai has odd many 
Z-paths from i because Definition [s] implies that /(ai) G s{i). If ai is already oddly connected to 5, 
then we are done and ai = cq. If ai is evenly connected to S C s(i)^ then we know that it must be 
oddly connected to 5 U {/(ai)} C s{i). There must exist another vertex /(a2) G s{i) \{S U {/(ai)}) 
to which ai is connected to for it to be evenly connected to s{i) as is required by Lemma [6j Because 
/(a2) G s{i) we know there exists a stepwise influencing path pi(/(a2)) (Lemma Isl) and we can 



extend that path by ai and /(ai) as was proven in Lemma 10 We move backward on this path 
and find the element a2. If a2 is oddly connected to 5, we are done and set cq = a2. Otherwise we 
can continue as was the case for ai until we find an element that is oddly connected to S. This 
element must exist since graph is finite and the Z corrections do not create any loops. We select 
Co — Cirri' Note that cannot be i because f{i) G 5 = s{i) H g{i) but /{am) ^ g{i). 

There is a Z-path from = cq to ai (we moved backwards along this path to find a^) and 
from ai to bn because of the way we selected ai. There also exists a Z-path from bn to every other 
bj such that < j < n, thus will also have a Z-path to every bj in 62, ... , bn}- Even more, 
because: 

dm > bn A bn > bo 

^ yjeA a^^f-\j) 

This completes the proof. □ 



Finally we could put together Lemmas^ and 14 



Lemma 15 (Equality of the penultimate SSF and optimal gflow layer). Let {G,I,0) be 
an open graph with flow [f^^f), corresponding SSF (5,^5) and optimal gflow {g,^g) such that 
\I\ = \0\. Then V{^' = V^' . 

Proof. Assume V-[^^ ^ V^^ we show how we can choose infinitely many different vertices from V{G). 
Due to Definition [g] we have \V{^^\ < \V-[^^\ and since V{^'' ^ V{^'' hence trivially V^' <t V^' and 
there must exist a vertex i in 

^i" \ ^\"- Then from Lemma [i] we have — O and using Lemma 



13 we obtain the following: 
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• gii) c O 



• f{i) e g{i) 

• g{i)ns{i) c s{i) 



and that there exists another vertex bo such that 



• 60 e Odd{g{i) n s(i)) 

• fipo)es{i)\g{i) 



These constraints together with an empty set A allow us to apply Lemma[T4| Lemma[T4|is constructed 
in such a way that whenever we can apply it to a (possibly empty) set A, it proves the existence of 
another set B such that that l^l < 1^4 U S| and Lemma ImI is applicable to the new set AVJ B. Thus 



it is possible to apply Lemma 14 infinitely many times and construct a subset of V{G) containing 
infinitely many vertices. This leads to a contradiction as G is a finite graph. □ 



2.3.2 Reducing the open graph 

The equality of penultimate layers of SSF and gflow might suggest that one could prove the equality 
of other layers simply by removing the last layer from the open graph and reapply the lemmas from 
the last section. However this would fail as the vertices in any layers can also use the output vertices 
in their correcting sets. Therefore we need to be careful which vertices we remove such that the 
reduced graph still have a gflow. 

Definition 9. // (G, /, O) is an open graph with flow (/, and corresponding SSF (5, ^5) then 
we call the open graph {G' ^I^O') a reduced open graph according to (5, -<s), where 

• R—{v^O\ f~^{v) G V^^} is the set of removed vertices. 

. G' ^ iy'.E') where 
V' ^V\R 
E' ^ E\{V X R) 

. o' ^ {v^^ yjO)\R 

We will omit "according to..." and call {G\I,0^) just reduced open graph when it is clear from 
the text which SSF is used for constructing it. An example of a reduced open graph is shown in 
Figure [6] 

As we saw in the previous section, we needed the fact that |/| = |0| to be able to prove that the 
penultimate layers of SSF and optimal gflow are equal. If we want to apply the same lemmas to the 
new reduced open graph, we need to guarantee that if we start with a graph where input size equals 
output size, the same holds for the reduced open graph. 

Lemma 16. Let {G'^I^O') he a reduced open graph of the open graph {G^I^O), then \0\ — |0'|. 

Proof. Let R be the set of vertices removed from G, then for every vertex i G V-^^ we have a 
corresponding unique vertex f{i) in R since Proposition [2] implies that s{i) C O and f{i) G s{i). On 
the other hand, for every vertex in R there exists a corresponding vertex in V-^^ from the definition 
of R. Therefore for every vertex v G R that we remove from O when constructing = (^1^^ UO)\R 
we add another vertex f~^{v) G V{^^ and it must hold that |0| = |0'|. □ 
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Figure 6: An example of an SSF reduced open graph (right) together with the original open graph 
(left). 

The next lemma is used later to construct a gflow of the reduced open graph from the gflow of 
the original open graph. 

Lemma 17. Let (G, /, O) be an open graph and A and B two sets in O such that Odd{B) H — 0. 
Then Odd{{A U S) \ (A n B)) n = Odd{A) n O^. 

Proof. There are altogether four different possibilities for a vertex v G to be connected to the 
sets A and B satisfying Odd{B) H = as shown in Figure u| 



1) 


V e Even{A) n Odd{A \B)^ve Odd{A n B) ^ 


V e Odd{B \A)^ 




V e Even{{A \B)U{B\ A)) 




2) 


V e Even{A) D Even{A \B) ^ v e Even{A n B) 


^ V e Even{B \ A) 




V e Even{{A \B)U{B\ A)) 




3) 


V e Odd{A) n Odd{A \B)^ve Even{A n B) ^ 


V e Even{B \A) ^ 




V e Odd{{A \B)U{B\ A)) 




4) 


V e Odd{A) n Even{A \B)^ve Odd{A n B) ^ 


V e Odd{B \A)^ 




V e Odd{{A \B)U{B\ A)) 





We see that every time v is evenly connected to A it is also evenly connected to {A \ B) U {B \ A) 
and every time v is oddly connected to A it is also oddly connected to {A \ B) L) (B \ A) . Because 
{A\B)U{B\A) ^ {AuB)\{Ar\ B) ans v is in it must hold that Odd{{A {JB)\{Af]B))r\0^ ^ 
Odd{A) n O^. □ 

We start by creating a function that will be proven to have the required properties of the gflow. 

Lemma 18 (Finding the reduced gflow function). Let {G,I,0) be an open graph with flow 
(/, -< f), SSF (s, -<s) and optimal gflow {g, ^g) such that \I\ = \0\. Let {G', I, O') be the SSF reduced 
open graph of {G, I, O) with the removed vertices set R, then there exists a function g' : 0"~" — > 

1. Vi G O'^ g'(i) n O'^ = g{i) n O'^ 

2. Vi G O'^ Odd{g'{i)) n O'^ = Odd{g{i)) n O'^ 
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Figure 7: The four possibilities for a vertex v G O to be connected connected to the sets A and B. 
The vertex v can be either oddly (a) or evenly (b) connected the sets A\B, B\A and AnB] oddly 
connected to A\B and evenly to B\A and An B (c); or evenly connected to A\B and oddly to 
5\ A and AnS (d). 



Proof. We start by noting that according to Lemma |3] we can create an optimal gflow ^g) of 
the open graph (G, /, O U V^') = (G, /, U i?) by restricting g to V{G) \ {V^' U V^') = O'^ and 
setting ^g—^g Wi^^ X ^o^^- construct our desired g' function from g. 

We consider i G O'^, if there exists a vertex j G RD g{i) then from the reduced open graph 
definition we have f~^{j) G V^^ . Also from Lemma 15 we have V^^ = V-[^^ and thus f~^{j) G y^^^. 
According to Proposition [2] this means that s{f~^{j)) C O. We have Odd{s{f~^{j)) D O'^ = since 
the only odd neighbours of s{f~^{j)) are either output vertices or the vertex f~^{j) G V^^ C 0^ 
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Now we define g^{i) = (^(z) U5(/ ^{j))) \ {g{i) ^s{f ^{j)))^ hence j g^ii)- Moreover Lemma 
implies that Odd{g'{{)) nO'^ = Odd{g{i))nO'^ . Also g\i)nO'^ = gli)^0'^ since s{f-^{j)) C 
Note that, since the new set will be constructed via a union of two sets we might add another vertex 
G i? to the set g\i). However, we can remove any such vertex k added to g\i) by applying the 
same procedure recursively. For every such vertex fc, it must hold that f~^{j) ^/ f~^{k) since 
k G s{f~^{j)) and Proposition [2] implies the existence of a Z-path from f~^{j) to f~^{k). Now 
we remove k via the above procedure i.e. defining g'{i) — {g'{i) U s{f~^{k))) \ {g\i) H s{f~^{k))). 
If this would add vertex j again to g^{i), hence there exists a Z-path from f~^{k) to f~^{j) and 
f~^{k) f~^{j) which contradicts the previous relation. This procedure will eventually terminate 
and remove all undesired vertices k ^ R since in the above procedure we never create any Z-path 
loops. □ 



We call a function which satisfies properties (1) and (2) of Lemma 18 the reduced gflow function 
of g. We can interpret these properties as saying that the gflow function g^ differs from the gflow 
function g only by the vertices in 0\ i.e. the other elements in the correcting set are left unchanged. 
As a gflow consists of a function and a partial order, we still need to define a valid partial order. 
The one that is most useful to us is such that it preserves as much relations as possible from the 
original gflow, hence the layering structures remain similar. 

Lemma 19 (Constructing the reduced gflow). Let (G, /, O) be an open graph with SSF (5, ^s)? 
gflow (g^^g) and {G'^I^O') a reduced open graph of {G' ^I^O'). If g' is a reduced gflow function of 
{g^^g), then (g^^g^) is a gflow of {G' ^I^O'), where 

yij eO'^ i^gj^i^g^j 

ViGO^^,VjGO^n^'(z) ^ i^g^j 



25 



Proof. We will show that {g\^gf) satisfies the three gfiow properties (Gl) - (G3) in Definition 
3| First property requires tha t if j G g^i)^ then i -<gf j. This is obviously true if j ^ O' . If 



j G g'{i) n O'^, from Lemma 18 we have j E g{i) which implies that i j because ^g) is a gfiow. 
Now according the definition of it must also hold, that i -<g' j. 

Now we consider the gfiow property (G2). For every j E Odd{g'{i)) it must be that j — i oi i ~<gi j. 
If j E then again this is obviously true because of the definition of ^^/. If j E Odd{g'{i)) H O'^ 
then we know that j E Odd{g{i)) and j — i oi i ^g j . According to the definition of ^^/, i -<g j 
implies that i j and we have that if j E Odd{g'{i)) then either i — j or i j. Thus the gfiow 
property (G2) is satisfied. 

Finally, we require for gfiow property (G3) that i E Odd{g'{i)) and as z E O'^ this is true because 
of the properties oi g' . □ 



We call the gfiow [g' ^ ^^/) from Lemma 19 the reduced gflow of (^, -<g). Similarly we can construct 
the SSF of the reduced open graph. Note that an SSF can only exist if the reduced open graph has 
fiow. Thus arises the need to prove the existence of a fiow on the reduced open graph, as is done in 
the next lemma. 

Lemma 20. // (G, /, O) is an open graph with flow (f^^f) and if {G' ^I^O') is the reduced open 
graph described in Definition^ then ^f), where 

• Vi G O'^ f'{i) = f{i) 

• ^//=X/ \[{V xR)U {V^^^ X O)] 
is a flow of{G',I,0'). 

Proof. It is sufficient to show that (/', -<fi) satisfies the fiow properties (Fl) - (F3) in Definition 
[2] and that /' is a function from O"^ to I*^. It is easy to see that /' : O"^' V \I. f acts by 
definition on O"-^ and 

VieO'^ /'(i) = /(i)ey\/. 

The graph has fewer vertices than G, therefore we need to show that all the vertices required 
according to the fiow function f are included in G\ i.e. ^i E O'^ it must hold that f\i) E V. 
According to Definition |9] every vertex v removed from the initial open graph (G, /, O) is chosen such 
that f~^{v) E 0^ Therefore it must be that every vertex j E such that f{j) ^ must be an 
output vertex in (G^ /, O'^). Because f\j) is not defined for outputs the vertices removed from the 
original graph G are not needed for f and f : O'^ V' \I. Hence we have that f{i) = f'{i) E 
for every vertex i E O^^ and : O^^ . 

Let R be the set of removed vertices as defined in Definition |9j The fiow property (Fl) states 
that i f'{i) and holds because: 

VzeO'^cO^ i^f{j)^f\j) ^ 

^ f'{i)iR A (z,/(z))E^/ ^ 

^ (z,/(z))E^/\y Xi?=^;, ^ 
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To prove that ^j/) satisfied flow property (F2) we need to show that for every j G if 
j G N{f\i)) then either j — i or i j. 

jeNifii)) jGiV(/(i)) 

^ j = i V i j ^ 

^ j = z V i j 
Finally the flow property (F3) i G N{f'{i)) holds almost trivially: 

z G O'^ C ^ z G 7V(/(z)) = A^(/(i)) 

□ 

Next we prove that the reduced gflow of an SSF is also an SFF. 

Lemma 21 (Constructing the reduced SSF). Let {G,I,0) be an open graph with flow (/, ^f) 
and SSF (5, ^5). // /, O^) is the reduced open graph, according to (5, ^s), then there exists an 
SSF {s\ ^5/) of {G\ /, O') such that {s\ ^5/) is the reduced gflow of (5, ^5). 

Proof Let R be the set of vertices removed from (G^I^O) to get {G'^I^O'). The reduced flow 



^fi) exists because of Lemma 20 Define (5', ^5/) to be the the SSF derived from this reduced 



flow. Assume (5', ^g/) is not a reduced gflow of (5, ^5), then one of the properties of Lemma 18 
should not hold, We show a contradiction in both cases. 
If the first property does not hold then 

3i G o'^ s.t. s\i^ n o'^ ^ s{i) n o'^ ^ 

3jGO^^n[(.(z)\.'(z))U(.^(z)\.(z))]^ 
{parity {a if = I A parity{Cf {f-\j))) = 0)V 
{parztyiCf {f-\j))) = I A parity{C^ {f-\j))) = 0) ^ 

parity{Cl{j)) parity {(f (j)) 

Hence by removing vertices and edges from the open graph (G, /, O) we must have changed Ci{f~^{j)) 
by an odd number to get Q {f~^{j))' 

We look at how removing the vertices in R from the open graph (G, /, O) changes Ciif~^U))- 
Removing a vertex v changes the number of Z-paths from i to f~^{j) if by removing it we also 
remove an edge in the Z-correction graph Gz- Let this removed edge be (A:,/), then Corollary [l] 
implies that / G N(f(k)) and v has to be either A:, / or f{k). Corollary [l] also implies that for v 
to have an outgoing edge in G^, f{v) has to be defined. Since / is not defined for output vertices 
and V ^ R C there cannot be any outgoing edges from v. Therefore v cannot be k as there is 
an edge from k to I in Gz- Also v cannot be / since again v cannot have an outgoing edge in G^, 
hence v would have to be the last element on the Z-path from i to f~^{j)^ which is f~^{j). This is 
impossible, as f~^{j) cannot be an output vertex. Therefore the only possibility is that v — f{k). 

Let V be the first vertex removed from G, such that Ci{f~^{j)) changes by an odd number. 
Hence all the paths from i to f~^{j) that disappear due to removal of v have to go through f~^{v). 
Therefore there must also exist an odd number of paths from f~^{v) to f~^{j). We know that 
because of Proposition |2| j G s{f~^{v)) and f~^{v) -<s j- On the other hand because of Definition [o] 
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it must also hold that f~^{v) G V{^% which together with Definition |4] implies that j G V^^ = O. 
This leads to a contradiction, because j has to be in O'^ C and cannot be in O. Therefore 
property (1) must be true for ^s')- 

Now we show that property (2) has to hold. According to Lemma [gI 



Odd{s'{i)) n o'^ = {i} = Odd{s{i)) n 

Because i G O'^ C O^, it must also hold that 

Odd{s\i)) n o'^ = {i} = Odd{s{i)) n o'^ 

and property (2) has to be true for {s' ^ ^s')- C 



In the specific case where |/| = |0|, it will follow from Lemma 21 that the unique SSF of a 
reduced open graph is the reduced gflow of the original SSF. 

Corollary 3. Let (G, /, O) he an open graph with an SSF (5, ^5) such that \I\ = \0\. If (G^ /, O') is 

the reduced open graph, according to (5, ^s)? of (G, /, O), then the unique SSF {s\ ^g/) of {G\ /, O') 
has the following properties: 

1. Vz G o'^ s\i) n o'^ = s{i) n o'^ 

2. Vi G o'^ Odd{s\i)) n o'^ = Odd{s{i)) n o'^ 



Proof. Because of Lemma [2Q| (G^ J, O^) has a flow. Since |/| = |0|, then according to Lemma 16 
|/| = \0'\ and hence (G'.I^ O') has a unique flow [3T]. Flow is required for the existence of an SSF 
according to Proposition [2} therefore there can exist only one SSF and because of Lemma 11_ this 



SSF has to satisfy properties (1) and (2). □ 
2.3.3 Moving back 

We have proven that the penultimate layers of SSF and optimal gflow are equal if |/| = \0\. Then 
we showed how to remove some vertices from the open graph and construct an SSF and optimal 
gflow on the new reduced graph. Both of them are reduced gflows, a property which we will use in 
this section to show that they preserve the layering of the gflows they were derived from. 

Lemma 22. Let {G,I,0) be an open graph with SSF (5, ^5) and gflow {g, ^g) such that {G\I,0^) 
is the reduced open graph of (G, /, O) with the removed vertices set R. For every reduced gflow 
{g^ -<g') of {G\ /, 00 such that V^' = V^' = O and V{^' = V^' it must hold that 

Vn>0 ULo n^^^ = U^iX^n (15) 



Proof We prove Lemma [22] by induction and showing first that Equation [TS] holds if n = 0, i.e. we 
need to prove that 

v,^'' = {v,^'uv,''n\R 

Lemma [2] tells that V^^' = and V^^ = O. Because the penultimate layers of SSF (5, ^5) and 
gfiow (g^^g) are equal we also have that V-[^^ — V^^ . Now we take the definition of O' from 
Definition [9] of the reduced open graph and substitute the appropriate sets: 

O' = iy^^ UO)\R ^ V^'^ = {V;^' U Vq^') \ R 
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Thus Equation 15 holds for n = 0. For the induction step we assume that Equation 15 holds for 
n = m — 1, i.e 

ur=o V = ur=on^n (16) 

and show that it holds for n = m. We use contradiction and assume that 

There are two possibilities: either 3i G V^^' \ V^^-^ or 3i G V^^-^ \ . We note that according to 

Lemma 
that 
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I j <^ i ^g' j if i G O'^ and j G O'^ . Because V^^' C O'^ for every m > we have 



^ 3j ev;^'' ng{i) s.t. i ^g, j ^ 

As can be seen above, both of the possible cases lead to a contradiction and hence it must hold that 

This completes the induction step and the proof itself. □ 

From the previous lemma we can construct a proof saying that every layer of a reduced gflow 
starting from the second to last one is equal to a layer of the original gflow. 

Corollary 4. Let (G, /, O) be an open graph with SSF (5, ^5) and gflow {g, ^g) such that {G\ /, O^) 
is the reduced open graph of (G, /, O) with the removed vertices set R. For every reduced gflow 
{g^ ^g') of {G\ /, 00 such that V^' = V^' = O and V{^' = V^' it must hold that 

Vn>0 Vn'' =V^_^^ 



Proof This follows trivially from Lemma [22j We have that if n > 0: 

UUV,"^' =UltlV,''^\R 

We can now subtract the elements of the second set from the first. 

Because Definition [o] of SSF reduced open graph we know that R C. O. Lemma [2] says that O TV^ 
hence we know that no element in R can be included in V^_^^ forn > and Vn ' = V^_^^. □ 

It turns out, that if the gfiow we have for the original open graph is the optimal one, then the 
reduced gfiow will be optimal for the reduced open graph. 

Lemma 23 (Constructing the optimal reduced gflow). Let {G,L,0) be an open graph with 
SSF (5, ^s) ci'^d optimal gflow (^, ^g). If (G\ /, O') is the reduced open graph of (G, /, O) then the 
reduced gflow ~<g') of (^, -<g) is the optimal gflow of (G\ /, O'). 
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Proof. First, because of Lemma [T9| {g\ g') has to be a gflow of {G\I,0^). Let us assume that 
-<g') is not the optimal gflow of (G^ /, O') and let (d, be the optimal one. Then according to 
Definition [6] 



3n > 0, i G O'^ s.t. i G \ T4 



A 



and since G U^^qT^^ 



from Lemma 22 



we obtain d{i) E U^^qT^^ ^ \ where the R 



is the set of vertices removed from the original graph. Now we know from Definition [6] that i is in 
V^_^'^ which according to Corollary |4 must be equal to Vn^' . This leads to a contradiction because 
i G \ Vn^' and thus {g\ ^^/) has to be the optimal gflow of {G\ /, O'). □ 



2.3.4 Proof of the optimality theorem 

We can now prove Theorem |2] by showing that the vertex layering of any SSF and an optimal gflow 
is exactly the same. Let (G, /, O) be an open graph with flow (/, ^/) such that |/| = |0|. Let 
(5, ^s) be the SSF obtained from (/, ^f) according to Proposition [2] and ^g) the optimal gflow of 
(G^I^O). According to Proposition [2] the last layer of any SSF is the set of output vertices. Lemma 
b^says that this is also true for the last layer of an optimal gflow, therefore V^^ = V^q^^ = O. The 
layers V-[^^ and V^^ are equal because of Lemma 



Now we need to show that layers and are equal for n > L We can construct a reduced 
open graph (Definition [9]) {G',I,0') from {GJ,0). We now consider the unique SSF ^5/) and 
reduced gflow {g' ^ -<gi) of {G\I,0^)^ which according to Lemma [23| is optimal. According to Lemma 
22j Vn^' = Ki+i every n > and because SSF is by Theorem [l] a gflow the same lemma also 
implies that V^^^ = Ki+i- Because of the way a reduced open graph is defined, we know that 
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|/| = \0^\ (Lemma 16). Thus we can again use Lemma 15 to say that = V^^' = V-[^^' = 



We can now take {G'^I^O') and find its reduced open graph to show using the same technique 
that V^' = This can be continued until we reach the empty layers, in which case we have 

considered all the layers according to (5, ^5) and (^, ^g). As every layer of (5, ^5) and (^, ^g) will 
be equal and (^, ^g) is the optimal gflow, the SSF of a flow of an open graph (G, /, O) is an optimal 
gflow if |/| = |0|, which proves Theorem [2j 



3 Compact circuits from signal shifted flow 

In the last section we presented an automated parallelisation technique for measurement patterns. 
However, when we translate those parallel measurement patterns back to the circuit model using 
the method described in [IT], we end up with quantum circuits with many ancilla qubits. More 
specifically, the new circuits will have the same number of qubits as there are vertices in the associated 
MBQC graph. Our next main result of the paper is a new scheme that explore the notion of circuit 
compactification introduced in [13j to remove all ancilla qubits introduced by the back-and-forth 
translation between the two models. We start by reviewing the notion of extended circuits, which 
is basically a re-interpretation of measurement patterns using circuit notation. Then we derive a 
set of rewrite procedures that combine the rewrite rules introduced in [13j in such a way that the 
SSF layering function (and, consequently, the optimised depth) does not change in the process of 
removing ancilla qubits from an SSF extended circuit. Finally, we introduce the algorithm that 
make use of rewrite procedures to completely rewrite an extended circuit until all ancilla qubits are 
removed. 
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Figure 8: Graph associated to the measurement pattern in Equation 17 and the corresponding 
extended circuit. 



3.1 Extended translation 

A straightforward translation method for measurement patterns, which we refer to as extended 
translation^ was introduced in [11]. This translation is inefficient, in the sense that it gives as many 
circuit wires as vertices in the original pattern (instead of inputs only). However its importance 
comes from the fact that its very easy to implement since the procedure to obtain an extended circuit 
is just a re-interpretation of the measurement pattern using the quantum circuit notation. Moreover, 
it will serve as a starting point to obtain more compact circuits for patterns with signal-shifted flow. 

Definition 10. Given a signal shifted measurement pattern with computational space (V^ /, O) and 
underlying geometry (G^I^O) with SSF (5, ^5). The corresponding extended circuit C with \I\ input 
qubits and \ /| ancilla quhits, is constructed in the following steps: 

1. Each vertex on the graph is translated as a circuit wire. 

2. The wires corresponding to are prepared in the sate |+). 

3. Each edge linking vertices i and j on the graph (command Eij) is translated as a CZij in the 
beginning of the circuit. 

4. Each dependent measurement [M^^'J^^^^ is translated as a gate J{—9i) in wire i followed by 
Controlled-X operator with qubit i as control and j G s{i) as the target. The layering respects 
the by replacing the control just after the J-gate and the target just before the J-gate of the 
next measurement command. 

5. Each correction on the output qubits (C — Pauli -X or -Z) is translated as a Controlled-C 
gates at with qubit i as control and j as target. The layering again respects ~<s by putting the 
control right after the J-gate and all the corresponding CX (introduced in Step 4) acting on 
qubit i. 

6. All the qubits in will be measured in the computational basis. 
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The obtained layering structure is referred to as £n, Jn, cind Cn, each containing only entangling 
gates, J -gates, and correction gates, respectively. We also divide slices Cn into Sn many smaller slices 
Cn^i, where Sn is the total number of J -gates in slice Jn- Each slice Cn^i contains all correction gates 
with control on qubit i s.t. Ji is in Jn- 

It is easy to verify that the above circuit implements the same operator as the measurement 
pattern (see also [H]). For clarity, in what follows we will refer to a CZij created in Step 3 above 
as Eij while keeping the notation CZij for those created in Step 5. Later on we will use the fact 



that, for a signal shifted pattern (Equation 14), a CXij or CZi-^ will be created in the corresponding 
extended circuit if and only if j G s{i) or /c G Odd{s{i))^ respectively Note that by construction, all 
gates associated to operators Eij are initially in slice £i (Step 3 in Definition 10), with £2^ -•-^£n all 



empty. However during the compactification procedure while we rewrite the circuit new gates will 
be added to these empty slices. Figure |8]-b shows the extended circuit of the following signal shifted 
measurement pattern with associated graph given in Figure |8]-a. 

3.2 Compactification procedures 

Compactification procedures can be described as a way of globally rewriting extended circuits in 
order to remove ancilla (non-input) wires. One way to achieve this is to rewrite the circuit to create 
J-blocks^ defined as follows. 

Definition 11. Consider a measurement pattern with computational space {V,I,0) and underlying 
geometry (G, /, O) with flow (/, ^j) and corresponding extended circuit C . We say there is a J-block 
in wires i and f(i) if the following set of conditions are satisfied (see Figure^a): 

1. The initial state of wire f(i) is |+). 

2. The gates sequence (E^f(^^^, JiiOi), appears in C . 

3. The only gate acting on the wire f{i) before CX^j(^) is E^fi^^y 

4. The only gates acting on wire i after E^f(^j^ are Ji{6i) and CX^j(^). 

5. After CX^j(^) gate, the qubit i is measured in the Z basis. 

Once a J-block is created (via circuit rewriting), one can use the identity in Figure [o] (J-gate 
identity) to remove one wire from the circuit. In general, extended circuits are not prepared for 
direct applications of the J-gate identity. In Definition 
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all corrections are translated as 



controlled-gates with control placed after J^ gate. Hence Condition 4 in Definition Vi is not satisfied 
in general. Moreover, since all E gates are initially placed in slice £1^ Condition 3 is not satisfied 
either for any open graph with more than two non-output qubits (see example in Figure |8]-b) . In 
order to create J-blocks in SSF extended circuits, we explore the relation between the E gates and 
the correcting gates C, since the latter are defined accordingly to the former through the stabilizer 
formalism. In other words, there is a direct relation between the gates in slice £1 and all other 
two-qubit gates in the rest of the extended circuit. In the case where we succeed in removing as 
many wires as there are non-output qubits in the graph, we say that the resulting circuit is in a 
compact form. We will refer to circuits in the compact form as compact circuits. 

Definition 12. Let C be the extended circuit of a measurement pattern with computational space 
{V, /, O). We say that C can be put into a compact form if there exists a sequence of circuit rewriting 
equations such that the J-gate identity (Figure^ can be applied \0^\ times. 
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Figure 9: This J -gate identity will be used repeatedly to simplify generic extended circuits. Note 
that the J-gate angles differ from each other by a minus sign. 
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Figure 10: Circuit identities in figures (a) to (c) are the rewrite rules from [13] (the circuit identity 
in Figure (b) is true only if qubit k is in the |+) state). The circuit identity in Figure (d) (Figure 
(e)) is obtained by multiplying CZij^ (CXi^) in both sides of the identity in Figure (a) (Figure (c)). 



A collection of circuit identities with the purpose of exploring the aforementioned relation between 



gates in extended circuits to create J-blocks was introduced in [13j (Figure 10). A compactification 
procedure for graphs with flow was provided and some simple examples of graphs with gflow explored. 
In this section we present a novel algorithm able to rewrite SSF extended circuits to put it into a 
compact form. In what follows we present the circuit identities that will be used in the algorithm, 
which we will refer to as Rewrite Procedures (RPs). We refer to the i wire of each RP as the target 
wire^ ji, as the correcting wires and finally k as the neighbour wire. Moreover, when we need 
to emphasize which RP we are referring to we also add a superscript to the wire label; for instance 
fc(^) indicates the neighbour wire of RP2. 



• Rewrite Procedure 1. The circuit identity in Figure 11 moves gates (E^j-^^ ...^ Ej^j^) past 
gates (CX^j^, CX^j^), adding m many gates CZij^ to slice C in the process. Using the 
rewrite rule in Figure 10 -d each of those E gates can be moved past gates (CX^j^, CXij^) 
in C, creating a new CZik each time the rule is applied (Figure [TT]-c). 



Rewrite Procedure 2. The circuit identity in Figure 
shoe £ with (CX 
use the rewrite rule in Figure 
gates {Ekj^,...,Ekj^_^) into (CX,,^, 



12 



replaces gates (Ej^j-^^ E^j 



m 



10 



^j^) in slice C, removing gate CXij^ in the process. We 
b for each pair {{Ej^kEkjJ, {Ej^-ikEkjJ} transforming 

n5 •••5 CXj^_ 



-Ijnj 



as depicted in Figure 12-b. The new CX 



gates can be pushed forward to the beginning of slice C, since it commutes trivially with 
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Figure 11: Rewrite Procedure 1. This RP is basically several applications of the rewrite rule in 
Figure [TOj-d. Although all E gates are drawn in the £ slice, it is sufficient that those E gates are 
placed just before the corresponding CX gates (that is, with no gate in between)). 



Using the rewrite rule in Figure 10 -e we can commute (CXj^j^, CX 



Jn — ljn . 



{CXij-^, . 

pre-existing CXij^ in C will cancel out, resulting in the circuit depicted in Figure 12 -c. 



past 

, CXij^_^) creating (n — 1) many new CXij^ in the process, which together with the 
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Figure 12: Rewrite Procedure 2. If Ls(A:) = Ls{i), slices {£,J,C) and {£' ^J' ^C) become the same; 
The rewrite procedure remains exactly the same. 
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Rewrite Procedure 3. The circuit identity in Figure 13 replaces gates {Ej^j^, Ej^j^) 
in slice £' with (CX^^ j^^t.), CXj^f(^]^^) in slice C. We use the rewrite rule in Figure 
b for each pair {{Ej^kEkf(k))^ ---^ {EjnkEkf(k))} transforming gates {Ej^k^ ---^ Ej^k) into 
^jif{k)^ --'^^^jrifik)) (Figure 13-b). The new CX gates can be pushed forward to the 
beginning of slice C, since it commutes trivially with CXj^f(^j^y Using the rewrite rule in 
e we can commute {CXj^f(^j^y CXj^f(^j^^) past (CX^j^, CX^j^), creating n many 
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CXifQ^-) in the process. Since n is even, all those CX gates will cancel, resulting in the circuit 
depicted in Figure 13-c. 



In order to apply the J-gate identity for the pair of wires (z, f{i)) of a SSF extended circuit, we 



need to rewrite it until all conditions in Definition [TT] are satisfied. In a SSF extended circuit, the first 
two conditions are trivially satisfied for any pair of wires (z, /(z)) and, therefore, we need to rewrite 
the circuit to satisfy the other conditions. To do so we analyse each qubit k in the neighbourhood 
of 5(z), classifying it according to three different cases: (i) Ls{k) < Ls{i)^ (ii) Ls{k) > Ls{i) and 
f(k) G s{i) and (iii) Ls{k) > Ls{i) and f{k) ^ s{i). This separation into cases is necessary for two 
reasons: First, the distinction between Ls{k) < Ls{i) and Ls{k) > Ls{i) is necessary because we are 
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Figure 13: Rewrite Procedure 3. If = Ls{i)^ slices {E^J^C) and {£' ^J' ^C) become the same; 

The rewrite procedure remains exactly the same. 



interested in keeping the J-gate parallelization introduced by signal shifting and hence we need a 
different procedure to deal with each case. Secondly, in the case where Ls{k) > Ls{i), we use the 



rewrite rule in Figure [10]- b which deletes E gates. Since condition 2 in Definition [TT] requires the 
existence of gates of form E]^f(^]^^^ we will treat differently cases where f{k) G s{i) and f{k) ^ s{i) to 
guarantee those E gates will not be removed from the circuit. As we show next, for each case one of 
the RPs can be applied if a set of prior conditions are satisfied. 

Proposition 3. Let i be a wire in a SSF extended circuit s.t. Ji is in some slice Jn- If there exists 
a wire k s.t. (i) Ls{k) < Ls{i) and (ii) the set of gates {E^j^, Ej^j^} (with ji, jm G s{i)) can 
be pushed to slice Jn^ then RPl ( Figure [Tip can be applied. 



Proof. Since Ls{k) < Ls{i)^ the gate belongs to a future slice Jm {m> n). Also, ji, jm ^ s{i) 
implies i {ji?---?jm} and hence the gates Jji,---,Jjm slices after Jji as well. Moreover, 

ji, •••5 jm ^ implies there exist operators in the measurement pattern, which are 



translated to the extended circuit as a CXij^, ...,CXij^ in slice Cn^i (according to Definition 10). 
Thus, if every gate Ej^j^ j E s(z), can be trivially pushed to slice we have exactly the scenario 
depicted in Figure [TI]-a. Therefore, the circuit identity in Figure IT can be applied. □ 



For reasons that will become clear in the section |3.3| where the algorithm to obtain compact 
circuits from SSF extended circuits is introduced, we need to consider a case which is slightly different 
from the scenario described in Proposition [sj In this case, condition (ii) in Proposition [s] is not 
satisfied because some of the E gates are in slice Cn^i but cannot be pushed trivially back to Jn- The 
only scenario where that could happen is if there exists CXi^k in slice Cn^i and some of the E^j gates 
are placed past it (and hence cannot be pushed trivially to Jn). In this scenario, as we show next, 
RPl can also be applied. 

Proposition 4. Let i be a wire in a SSF extended circuit s.t. Ji is in some slice Jn- If there exists a 
wire k s.t. (i) k G s(i) and (ii) gates {Ej^j-^^ ...^ Ej^j^} (with ji^ ...^ jm G s(i)) can either be all pushed 
to slice Jn or just some can be pushed to slice Jn and the other E gates are placed in Cn^i just after 
CXij^, then RPl ( Figure [77p can be applied. 



Proof. The proof of this proposition is trivial due to its similarity to Proposition [3j First, note that 
k ^ s{i) implies Ls{k) < Ls{i) and hence condition (i) in this proposition is equivalent to the one in 
Proposition |3j Therefore, if all gates in the set {Ej^j-^^ ...^Ej^j^} can be pushed to slice Jn we have 
exactly the conditions in Proposition [3| and there is nothing to prove. The other possibility is when 
a subset of the set {Ej^j-^^ Ej^j^} can be pushed to slice Jn but gates in the complementary subset 
are placed in slice Cn^i, after gate CXik (which exists since k ^ s{i)). It is easy to note that it will 
not prevent the application of RPl, since the only gate in Cn^i that could be placed in between the 
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Figure 14: Removing even many CZij from a SSF extended circuit. See Lemma 24 for more 
information. 



E and CX gates used in RPl, namely CXik, is assumed in condition (ii) to be placed before the 
gates of the form E^j. □ 

In the next Lemma we show the interesting effect of applying RPl to SSF extended circuits. 

Lemma 24. The application of RPl to a pair of target and neighbour wires (i, k) of a SSF extended 
circuit removes all CZik from the circuit. 



Proof. According to Definition 10, all CZ gates with control in wire i are placed in slice Cn^i and 



hence we only need to show that all CZij^ in that slice are removed. Let us divide the analysis into 
two cases: (i) k G Odd{s{i)) and (ii) k ^ Odd{s{i)). Consider the first case. Since k G Odd{s{i)), 



there exists a Zf' in the measurement pattern and hence there exists a gate CZij^ in slice 



For 



this case, the index m in RPl is an odd number and therefore the application of RPl creates odd 
many CZi^ gates in slice Cn^i- Note that those CZij^ gates can be created in different parts of Cn^i 
(depending whether Proposition [s] or [4] is satisfied) such that they cannot be grouped together 
trivially (like in Figure [Mj-a). It is easy to verify that the only possible scenario for that to happen 
is when there exists CXik in Cn^i in between the created CZik gates. However, those CZik gates can 
be grouped together by moving all CZik to one side of the troublesome CX^k using the identity 

CZikCXik\tl^)i\(l))k = CXikCZik{Zi\il;)i)\(l))k 

where and are arbitrary quantum sates and Z is the Pauli operator (see Figure [T4|-b) . This 
way, since there exist even many CZi^ in sequence and even many CZ gates equal the identity, 
all CZik can be simply removed from the circuit. The created single-qubit gate Zi can be pushed 
forward to the end of the i wire, since it commutes with all gates in Cn^i- Moreover, since the final 
measurement is onto the Z-basis, the operator Zi has no effect in the measurement statistics and 
therefore can be removed from the circuit without changing the computation being implemented by 
the circuit (Figure 14 -c). An equivalent analysis applies for the second case, where k ^ Odd{s{i)). 
For this case RPl create even many CZik in slice Cn^i and there is no pre-existing CZik in that slice. 
The same identity can be used to group together all created CZik-, which cancel out since there are 
even many of those. Since both created and pre-existing CZik are removed from the circuit by the 
application of RPl, the Lemma holds. 

□ 

Now we analyze the cases where Ls{k) > Ls{i)^ that is, with Jk gate blocking all gates of form 



Ekj to be pushed past the correction gates (necessary to satisfy condition 3 in Definition 11). Since 
in those cases the corresponding RP will delete E gates and we do not want to delete E gates of the 
form Eiji^i) (see Definition 11), we analyze separately the case where f{k) G s{i) (Proposition [5]) 
and the case where f{k) ^ s{i) (Proposition [g]). 
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Proposition 5. Let i be a wire in a SSF extended circuit s.t. Ji is in slice Jn- If there exists a 
wire k such that the following conditions are satisfied: (i) Ls{k) > Ls{i); (ii) f{k) G s{i); (Hi) gates 
Ekj^,...,Ekj^ (with G s{i)) can he trivially pushed to some slice Jn' (containing Jk); and 



(iv) the E]^f(^]^^ gate is the first gate acting on wire f{k), then RP2 (Figure 12) can be applied. 



Proof Since Ls{k) > Ls{i), n' < n. In slice Cn' there is a CXj^f(^j^y, since f{k) G s{i), we define 
jn — /(^)- Also, ji, ^ s{i) implies that i ^5 {ji, jVi} and hence the gates Jj^^ ...^ Jj^ are in 
some future slices (compared to n). Moreover, it also implies the existence of the operators ...,Xj^ 
in the measurement pattern, which are translated to the extended circuit as CXij^^ CXij^ in slice 
Cn^i (containing all the CX and CZ with control on qubit i). On the other hand, Ls{k) > Ls{i) 
implies both k ^ Odd[s{i)) and k ^ s(i) and therefore there is no CZij^ or CXij^ in the circuit. Thus, 
if every gate E^j^^ ...,Ekj^ can be trivially pushed to slice Jn' and £'/c/(fc) is the first gate acting 



on wire j^, we have exactly the scenario depicted in Figure 12 -a. Therefore, the circuit identity in 



Figure [12] can be applied. □ 

Proposition 6. Let i be a wire in a SSF extended circuit s.t. Ji is in slice Jn- If there exists a 
wire k such that the following conditions are satisfied: (i) Ls{k) > Ls{i); (ii) f{k) ^ s{i); (Hi) gates 
Ej^j-^^ ...^Ekj^ (with ji, ^ ^{^)) can be trivially pushed to slice Jn' (containing Jk); and (iv) the 



Ekf{k) gate is the first gate acting on wire f{k), then RP3 (Figure 13) can be applied. 



Proof. Since Ls{k) > Ls{i)^ then < n. In slice Cn' there is a CX^fi^^^ with f{k) ^ s{i). Also, 
jii •••ijn ^ s{i) implies that i {jii •••ijn} and hence the gates Jj^, Jj^ are also in some future 
slices (compared to n). Moreover, it also implies the existence of the operators ...^Xj^ in the 
measurement pattern, which are translated to the extended circuit as CX^j^, CXij^ in slice Cn^i 
(containing all the CX and CZ with control on qubit i). On the other hand, Ls(k) > Ls{i) implies 
both k ^ Odd{s{i)) and k ^ s{i) and therefore there is no CZi^ or CXi^ in the circuit. Thus, if 
every gate E^j-^^ E^j^ and £^A:/(/c) can be trivially pushed to slice Jn' and £'/c/(fc) is the first gate 



acting on wire /(A:), we have exactly the scenario depicted in Figure 13 -a. Therefore, the circuit 



identity in Figure [TS] can be applied. □ 

In what follows we analyse some properties regarding the interplay between the RPs that will be 
crucial for the compactification algorithm for SSF extended circuits. 

Lemma 25. Let C be an extended circuit obtained from a signal shifted measurement pattern. Then, 
if a rewrite procedure removes a gate Ej^ from the circuit, it will not be required for the application 
of any other rewrite procedure. 

Proof. Note that only RP2 and RPS can delete a Ej^ gate from the circuit. The proof is divided 
into three parts: (i) we show that the Ej^ is never a fiow edge and, therefore, it will not be required 
for the application of a J-gate identity; (ii) it will not be needed for the application of RPl and (iii) 
it will not be needed to create a CX (by using RP2 or RPS) in the next layers. Suppose Ej^ is a 
fiow edge such that k = /(j), with j being any of qubits (ji, j^_i) in RP2 or (ji, ...,jn) in RPS. 
The other way around, that is, j = /(fc), is not possible since f{k) is separately identified in both 
RP2 and RPS. Since f{j) G s{j) for all j G O^, we have j ^5 k. By construction, RP2 and RPS are 
applied if and only if k i- Moreover, since j G s{i) in both RP2 and RPS, it holds that i -<s j. 
Putting everything together we have i ^5 j ^5 k ^5 i, which is a contradiction. Therefore, Ejk is 
never a fiow edge. 

In this second part of the proof, we show that if Ej^ is deleted by RP2 or RPS, it will not be 
required for the application of RPl in any future step. Using the notation defined earlier, the fact 
that RPl is applied after RP2 or RPS can be written as i^^^ ^5 i^^^ or i^^^ ^5 i^^\ where i is the 
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target wire of the corresponding RP. This relation will be used in the rest of the proof. We show that 
and fc^^^ cannot be any of wires ...^jn^ (denoted simply by in the rest of this proof) or 
k^^^ and hence the gates E-{^2)^{2) or Ej(3)f.(3) cannot be the same as £"^(1)^(1), which is what we want 
to prove. Let us show for k^'^\ Since fc^^^ ^5 z^^^ and i^-^^ -<s k^'^'^ (true for both Propositions [s] and [ij 
which state the scenarios where RPl is applied) we have A:^^^ ^5 z^^^ ^5 i^^"^ k^^"^ and therefore 
k^'^^ cannot be the same wire as k^^\ Now assume wire j(^) is k^'^\ In RP2 we have A:^^) z(2) (by 
construction) and in RPl it holds that i^^'' ^5 since the target of a correcting gate is always 



placed before the J-gate acting on that same wire (Definition 10). Putting everything together gives 

since by assumption j^^"^ is k^'^\ it gives fc^^^ -<s fc^^^ which is a 
contradiction. The proof is the same for k^^\ Therefore, if Ejk is "consumed" by RP2 or RP3, it 
will not be required for the application of RPl in any future step. 

Finally for the third part note that an E gate will be removed to create a gate only through 



identity in Figure [TQ]-b. The pair of such E gates is always associated with a flow and a non-flow 
edge when used in RP2 or RP3. Therefore, there are two possible pairs of E gates using the non-flow 
edge Ejk'. EjkEj^f(^j^^ (which would create CXjfQ^^) and EjkEjf(^j^ (which would create CXj^f(^j^). 
We want to show that whenever we "consume" Ejj^, it will not be required to create a different CX 
than the one already created. Suppose we use Exy to create CXr^f^y). This is the case when y i 
and z X, if X is a correcting wire and y is a neighbour wire in either RP2 or RP3. Similarly, if x is 
the neighbour wire and y is the correcting one then x ^5 z and y i (easily obtained by relabelling 
the wires in RP2 and RP3). Since the first pair of conditions is inconsistent with the last pair of 
conditions (j ^5 i and /c i), we conclude that once a given E is used to create a CX, that E 
would not be required to create a different CX in any further step of the algorithm since the partial 
order ^5 is not changed by the application of any RP. □ 

Corollary 5. Let CXij he a gate in an extended circuit created by the application of a rewrite 
procedure. The E gate '^consumed^^ to create this CXij can he univocally determined. 



Proof. It follows from the third part of the proof of Lemma 25 that if a gate of the form Ejk is 
deleted in RP2 or RP3 we will have one of the following cases: the gates EjkEj^f(^j^^^ are used to create 
CXjf(^]^^ or the pair Ejj^Ejf(^j^ are used to create CXj^f(^jy The CX that is created by "consuming" 
the Ej}^ depends exclusively on the relation between vertices j, k and f{j) in respect to the partial 
order ^5. Since the partial order ^5 is not changed by the application of any RP, there is a one-to-one 
correspondence between the deleted E gate and the newly created CX gate. □ 

The process of choosing the correct RP to be applied is summarised in Algorithm [2j which will be 
used as a subroutine in the compactification algorithm for SSF extended circuits (Algorithm [s]) . The 
modified conditions in Lines 5 and 7 are based on the fact that if a CZ required for the application 



of RP2 or RP3 is not in the circuit, the CX that it would create will be (see Lemma 25), allowing 
the application of future RPs. In all RPs we assumed that some E gates could be moved to the 
beginning of a given slice J7 in a trivial way. It will not be true in general for SSF extended circuits, 
since there might exist several other gates in the extended circuit such that the aforementioned E 
gates could not be moved trivially to J . In other words, the conditions in Propositions [3] to [6] would 
not be satisfied. The algorithm in the next section address exactly this problem: it provides an 
ordering where Algorithm [2] never aborts. This ordering is the global structure coming into play, 
since it is related to SSF and flow of the graph. 
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Algorithm 2: This algorithm decides which Rewrite Procedure must be apphed by analysing 
how a given qubit k is connected with the correcting set of another qubit i 
Input: Wire labels i and k. 

Output: A description of which Rewrite Procedure must be applied, 
begin 

Read i and k. 

if all conditions in Proposition^ or Proposition^ are satisfied then 
I Apply RPl; 

if all conditions in Proposition^ are satisfied or conditions (Hi) and/or (iv) in Proposition 

are not satisfied but there exist a sequence of {CXj^j^^ CXj^_^j^) gates that can be 
pushed trivially to Cn^i (slice containing all the CX and CZ with control on qubit I) then 
I Apply RP2; 

if all conditions in Proposition^ are satisfied or conditions (Hi) and/or (iv) in Proposition 
1^ are not satisfied but there exist a sequence (CX^^j^j^.), CXj^f(^]^^) gates that can be 
pushed trivially to Cn^i (slice containing all the CX and CZ with control on qubit i) then 
L Apply RP3; 

otherwise Abort 



3.3 Obtaining optimised compact circuits 

In this section we explain how the compactification algorithm for SSF extended circuits (Algorithm 
^ works and show some examples. The goal of Algorithm^ is to create \0^\ many J-blocks (see 



Definition 11) in a SSF extended circuit and then apply the J-gate identity for all J-block, removing 
\0^\ many wires from the extended circuit. Therefore, the output of Algorithm [s] is the compact 
form of the inputed SSF extended circuit. One of the main differences between SSF and other gflows 



is the notions of stepwise influencing path, introduced in Section |2.2[ The next lemma use some 
properties of stepwise influencing path to relate the SSF correcting set to the partial order of flow. 
This relation will play an important role in Algorithm [3] since it gives an appropriate ordering for 
the application of the RPs. 

Lemma 26. Let {G,I,0) be an open graph with flow {f,^f) and SSF (5,^5). Then, for all 
V G N(j) \ {i}, where j G s(i), it holds that v yf i. 

Proof First suppose j = /(z). Then, by flow definition (Definition [2]), v yf i. Now suppose j ^ /(z) 
then from Definition |2j v >-/ f~^{j)- By Lemma [sj there exists a step- wise influencing path passing 
through f{i) and f {j)^ namely pi{j)- It follows from Corollary [2j that there exists f~^{k) such 
that k G s{i) and f~^{j) G N(k) and, consequently, f~^{j) ^/ f (k). Lemma 2.14 allows us to 
repeat this process to find the previous vertices in the path pi{j) until we reach vertex i. Hence we 
conclude that f~^{j) i and therefore v f~^{j) ^- D 

Before running Algorithm |3j the slices Cn^i must be arranged in the extended circuit from right 
to left respecting the order imposed by (see for instance Figure [s]). This can always be done 
since the set of gates in a given slice Cn^i commutes with the gates in an other slice c^j, for any 
valid j. The algorithm starts with a for loop that runs for one SSF layer at a time, starting with 
the input layer and moving onwards until the last layer of non-output qubit is considered. At each 
iteration n of the for loop, a rewrite procedure is applied to each pair of qubits (i, k) s.t. Ji is in 
J^n and k G N{s{i)). The two foreach loops and the two while loops define the order in which 
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Algorithm [2] will be called. This ordering, which is the inverse of the order given by ^j, assures that 
Algorithm ]2| will never abort when it is called by Algorithm [3| 



Algorithm 3: Transforms a SSF extended circuit into a compact circuit (Definition 



Input: SSF extended circuit. 

Output: Compact form of the SSF extended circuit, 
begin 

D ^ m^yii^G{Ls{i)} 
for n — 1 to D do 



Let Wn {an^i 



5 •••5 ^n,mn 



} be the set of the wires with J-gate in slice Jn 



maxievy„{^/(«)} 
Let Sl^i^ ^ miniewALfii)} 
while S^^^ > Sl^i^ do 

foreach an,i G Wn such that Lf{an,', 

Let NSTEP ^ maxfcg{^(,(„„^))\{„^^}}{L/(fc)} 
while NSTEP > S"^^^ do 

// That is, for all k E {N {s{an,i)) \ {o,n,i}} (Due to Lemma 
foreach k G {N{s{an,i)) \ {an,i}} such that Lf{k) = NSTEP do" 

Run Algorithm [2] 
NSTEP ^ {NSTEP - 1) 



26) 



{S„ 



I) 



Remove from the circuit all gates CXij placed in any slice Cn such that j ^ f{i) 
Apply the J-gate identity for all pairs of wires {z, f{i)} such that i G . 



Lemma 27. Algorithm^ never aborts when called by Algorithm\^ 

Proof. Let us start by showing that the algorithm works for the first iteration of the for loop {i.e., 
n = 1) and then we explain why it will work for all other layers. In what follows we use the notation 
Qn^i to represent a target wire i in layer n; note that there might exist more than one such wire in 
the same layer. Let S^^^ i'^min) maximum (minimum) value of Lf{an^i) for an^i G Wn (as 

defined in Lines 5 and 6 in the algorithm), where Wn is the set of the wires with J-gate in slice 



Jn. According to Lemma 26, for any ai^i such that Lf{ai^i) — ^niax^ Qubit k connected 

to a qubit in s{ai^i) it holds that k ^ Wi. In this case, condition (i) in Proposition [s] is satisfied. 
Since we are in the first SSF layer, all E gates required for the application of any of those RPs can 
be trivially pushed forward to slice J7i, hence all conditions in Proposition [3] are initially satisfied. 
After the application of one or more RPl, it might happen that some E gates get stuck in between 
two CX gates in slice ci^i (which happens when a neighbour wire k is itself in both N(s{ai^i)) and 
s{ai^i) sets). In this case the conditions in Proposition |4] are the ones satisfied. In both cases, only 
RPl can be selected by Algorithm [2| Therefore, Algorithm [2] would not abort for all target qubits 
ai^i such that Lf{ai^i) — Sj^ax- ^ consequence, since only RPl is applied, all gate Ekp^ such that 
P ^ s{ai^i), are moved past the correction slice ci^^, for all ai^i such that Lf{ai^i) — S^ax- 

Now consider qubits ai^i such that L f{ai^i) — Sj^ax ~ ^^at is, the second iteration of the first 



while loop. For k G N{s{ai^i)), Lemma 26 implies that either k ^Wioi {{k G Wi)A[Lf{k) = S^^J} 



is true. For k ^ Wi, the procedure to be applied is RPl, similarly and by the same reasons as 
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in the previous iteration. For the other case the apphcable rewrite procedures are either RP2 or 
RP3, since in Ji imphes Lf{k) — Lf{ai^i). Note that in the previous iteration (qubits ai^i s.t. 
Lf{ai^i) = S^ax^ lil^^ ^he wire k being analysed now), all Ef(^]^^y^ for any v ^ k were moved past the 
corresponding correction slice and hence Ef(^j^^j^ is the first gate acting on the f{k) wire (as required 
by both Propositions [5] and [g]) . Moreover, when the circuit identity in Figure [Toj-b creates the new 
CX gates (first step in both RP2 and RP3), it can be pushed forward to ci^i because the only gate 
in between is CXp.f(^j^^^ which commutes with the CX gates we want to push forward. Therefore, 
since all conditions in either Proposition [5] or Proposition [6] would be satisfied for all qubits ai^i such 
that Lf{ai^i) — S^ax ~ ^^e corresponding RP can be successfully applied. 

To complete the proof by induction for the first iteration of the for loop, assume Algorithm [2] 
has not aborted in the first q iteration of the first while loop. Then in the (g+ 1)^^ iteration, qubits 
ai^i such that Lf{ai^i) = {S^ax - q) ^ ^min^ where S}^-^ = min^i .Gi^i{^/(tti,i)}, are considered. 
If a given qubit k is neighbour of a qubit in s{ai^i)^ then Lemma [26| implies that either k ^ Wi 
or {{k G Wi) A [(S'^aa; ~ Qi) < Lf{k)} is true. Here the same analysis as before applies; First, 
if ^ VFi, Algorithm [2] will apply RPl, moving the corresponding E gates past slice ci^i. On 
the other hand, if {{k G Wi) A [{S^^^ — q) < Lf{k)}, either the conditions in Proposition [s] or 
Proposition |6] will be satisfied, since all E gates acting on qubits in the sets s{ G Wi 

s.t. Lf{ai^y) > {S^^^ — g), were moved past slice ci^y in previous iterations (with the obvious 
exception of gates of the form The procedure continues until there is no qubit left in Wi to 

be considered. Therefore, Algorithm [2] never aborts when called during the first iteration (n = 1). 

Now let us analyse how the E gates placement changed in the circuit after the n^^ iteration of 
the for loop where Algorithm [2] applied a RP each time it was called by Algorithm [3j First, note 
that all E gates acting on qubits in 5(z), for i G Wn^ were moved to £n+i (via RPl) or transformed 
into a CX (via RP2 or RP3) and then moved to the same slice, with the exception of gates of the 
form which remains in £n and will later be used for the application of the J-gate identity. 

Moreover, Lemma [25] guarantees that if an E gate is deleted in a given iteration of the algorithm, it 
will not be required in any future step of the algorithm. 

Note also that those E and CX gates were moved to the £n+i slice in a specific order, namely 
the inverse of the order induced by Since this ordering is a property of the associated graph, it 
does not change during the run of the algorithm. It means that for any given iteration n of the for 
loop, the ordering of the two-qubit gates required for the application of a RP will be in agreement 
with the order Algorithm [2] will be called by Algorithm [sj Hence the condition (required for the 
application of any RP) that the required two-qubit gates can be pushed trivially to a given slice is 
always satisfied. Therefore, at each iteration of the for loop, it arranges all two-qubit gates necessary 
for the next iteration in the same order it will be called. 

To complete the proof by induction, let us assume that in the first {p — 1) iterations of the for 
loop Algorithm ^ has not aborted, and then we show that it will not abort in the p^^ iteration. Let 



us analyse the first iteration of the first while loop for n — p. According to Lemma 26, for any a^^i 
such that Lf{ap^i) = Smaxi if a given qubit k G N{s{ap^i)) then either k G Wm (if there exists m < p 
s.t. S^^^^ > Smax) or the Jk gate belongs to a layer JT^/ s.t. > p. If the latter condition is the 
case, then by definition only RPl is applied to ap^i. On the other hand, if A: G Wm, RP2 or RP3 
must be applied. Since by assumption Algorithm [2] has not aborted in any previous iteration, the 
conditions for the application of the aforementioned RPs are satisfied and Algorithm 2 will not abort 
in the current iteration. 

Now assume Algorithm [2] has not aborted in the first q iteration of the first while loop for 
n — p. In the [q + 1)^^ iteration of this loop, qubits ap^i such that Lf{ap^i) — {Smax — q) ^ ^min^ 
where S^- = miiia^ .^Wp{L f{ap^i)} ^ are considered. For a given k neighbour of a qubit in s{ap^i), 
Lemma 26 implies that one of the following statement holds: (i) k G Wm A [S^^^^ > {Smax — q)] 
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for some m < (ii) {k G Wp) A [{S^ax — q) < Lf{k)] or (iii) k G W^'^ for > p. Here the same 
analysis as before applies; For cases (i) and (ii) RP2 or RP3 is applied and for (iii) RPl will move 
the corresponding E gates past slice Cp^i. Since the same order is respected throughout the algorithm 
(the inverse of the order defined by flow), the conditions for the application of the aforementioned 
RPs are satisfied and Algorithm 2 will not abort in the current iteration. This procedure is repeated 
until there is no more elements in the set Wp, which concludes our proof by induction and proves 



the Lemma. 



□ 



Theorem 3. Algorithm\^ outputs a compact circuit. 

Proof. We show that Algorithm [s] creates \0^\ many J-gate blocks (Definition 12) and then removes 
\0^\ many wires from the circuit (using the J-gate identity), yielding a compact circuit. The 
first two conditions in Definition [TT] are trivially satisfied for all SSF extended circuit. In Lemma 
27 we proved that Algorithm [s] moves to a future slice (or removes) all gates E acting on qubits 
j G s{i) (Vz G O^) except gate and hence condition 3 in Definition is also satisfied. Now 

we prove that condition 4 is satisfied. Note that all non-input qubit are initialised in state |+). 
According to SSF construction (Proposition [2]), for all qubit j G s{i), there exists f~^{j). By 
definition / : and therefore every qubits in 5(z), for all i G O^, starts in the |+) state. 

Therefore, all remaining CXij such that j ^ /(z) can me removed from the circuit without changing 
the computation, since CXij\-\-)j — (remember that the only gate E acting on wires j G s{i) 
left behind by the algorithm is E'^j^^)). Therefore, after the step in Line 17, condition 4 in Definition 
12 is also satisfied. Since all non-output wires in an extended circuit are measured in the Z basis by 
construction, condition 5 is also satisfied and hence Algorithm [s] has created \0^\ many J-blocks 
after the step in Line 17. Finally, in Line 18 the J-gate identity is applied to all J-blocks, resulting 
in a compact circuit. □ 

We complete this section with two examples using the above algorithm. In the first one we start 
with a generic quantum circuit, translate it to the MBQC model, find the SSF extended circuit of 
it and then apply Algorithm [3] to put the circuit in the compact form. In the second example, we 
analyse an example with Pauli measurements and show that if we want to keep the parallelization 
introduced by it we can no longer find the compact form of the extended circuit. 



Example 1 



In this first example we apply our method to optimise the circuit in Figure 16 -a. Using the method 



from |llj, we translate the circuit in Figure 16 -a to the following measurement pattern, resulting in 



the graph in Figure 15 a. The difference in depth between the flow and SSF is shown in Table 18 



Method 


depth 


partial order 


Flow 


5 




SSF 


2 


1,4,7^, 2,5 



Xg^e M^' Xg"^ z'^" zi" x^" x^" xi" Zg'i x*i Mp mI^ mI^ 

-E'78 -E'67-E'57-E'56 -E'45 -E'37 -E'35 -E'25 -E'24-E'23 -E'12 



(18) 



The extended translation (Definition [TO]) of the measurement pattern in Equation 18 gives the 
circuit in Figure 16-b. Let us apply Algorithm [3] to find the compact form of that extended circuit. 
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7 8 

(a) 



Figure 15: Open graph associated to the example in Section 3^ 



We start with the first SSF layer of circuit in Figure 16-a. For n = 1 we have Wi — {1,4,7}, 
5^^^ = L/(7) = 5 and 5^-^ = L/(l) = 1. Thus, we start with qubit 7 with A^(5(7))\{7} = {0} and 
therefore no rewrite procedure is applied. The value of 8}^^^^ is decreased to 2 (there is no qubit i s.t. 
Lf{i) = 3 or 4). Now we analyse qubit 4 because L/(4) = 2 with 7V(5(4))\{4} = {2, 3, 5, 6, 7}. For 
qubits in this set it holds that 2^j5^/7^/3,6. The maximal value of Lj is stored in NSTEP. 
For both qubits 3 and 6, Algorithm [3] applies RPl, moving CZq^ and past CX45. 

Now NSTEP ^ NSTEP - 1 and the qubit considered is qubit 7. For qubit 7 the algorithm 
applies RP2, transforming CZer, and into CXes, CX^^ and CX38 respectively and 

moving those new CXs past the correction structure, removing CX48 from the circuit. Now consider 
qubit 5 and move CZ^q and CZ53 past the correction structure, using RP2. The same rewrite 
procedure is applied for qubit 2, moving CZ23 and past the correction structure. The rewritten 

circuit is depicted in Figure [T6]-c. 

In the next step, the value of S}^^^ is decreased to 1, where the qubit to be considered is qubit 1. 
7V(5(1))\{1} = {2,3,4,5,6,7}. For qubits in this set it holds that 4^/2^/5^/7^/ 3, 6. For 
6, Proposition [3] is satisfied and hence Algorithm |3] applies RPl, moving CZ^^ past CX15. Then 
for A: = 3, Algorithm [s] applies RPl (Proposition [s]), moving £"35 and £^23 past CX12 and CX15. 
Next, RP2 is apphed to = 7, moving CX^^ and CX^^ past CX13 and CX15. For = 5, we have 
that k G s{i) (condition (i) in Proposition [4]) and £35 in slice ci^i after CX15 gate ((condition (ii) in 



Proposition [4]). Therefore, for this case RPl moves £35 and £25 past CX12 and CX13, creating two 
(7Zi5 separated by CX15. By Lemma 24, we know those two CZ15 can simply be removed from the 
circuit without changing the computation being implemented. For k — 2^ Proposition [3] is satisfied 
and hence RPl moves £25 and £23 past CX13 and CX15. Finally, for qubit 4, RP3 transforms £24 
into CX25 and move it past CX12 and CX15, removing the latter in the process, resulting in Figure 
[Md. 

In the second SSF layer we have n = 2, 1^2 = {2, 5}, 5^^^ = L/(2) = 3 and 5^-^ = L/(5) = 4. 
Thus, we start with qubit 5. Since A^(5(5))\{5} = {7}, Algorithm [s] applies RP2, moving CXq^ past 
CX56 and CX58 while removing CX58 in the process. After that, the value of NSTEP decreases 
and qubit 2 is considered. A^(5(2))\{2} = {5,7} and since £/(5) < Lf{7)^ qubit 7 is considered first. 
In this case Algorithm [3] applies RP3 moving gates CX38 and CXqs past the correction structure. 
Now for qubit 5, RP2 is applied transforming CZ53 into CX^e- The created CX is moved forward 
past the correction structure, removing CX26 from the circuit. The resulting circuit is shown in 



Figure 16 -f. The command in Line 17 of Algorithm [3] removes the remanning undesired gates using 
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Figure 16: A complete example of global circuit optimization. In each figure the shaded (blue online) 



gates are the new gates obtained through the RPs (see Section 18 for more information. 
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the identity trivial CXij\-\-)j = The resulting circuit is depicted in Figure 16 -g. Finally, the 

J-gate identity is applied for every qubit i G resulting in the optimised compact circuit shown in 
Figure fTBlh. 



Example 2 

In this example we explore what is the role of Pauli measurements in the context of compactification 
procedures. Let us start by reviewing how the correction operators are modified when applied to a 
qubit measured with and 7r/2 angles (corresponding to Pauli measurements): 

M^Xt = Ufzi (19) 
M^Xt = mO (20) 

Note that in the MBQC context, both substitutions clearly might reduce the depth of the computation: 



Equation 19 substitute Zf for Xf^ which can then be removed using signal-shifting, and Equation 
20 simply deletes the existing dependency. However, in the associated extended circuit, these 
optimisations can not be implemented alongside the compactification procedure developed in this 
section (Algorithm [s]) , forcing one to choose between optimization in time or memory. To see why. 



note that both Equations [19] and [20] change the correcting structure for every qubit encoded in the 
index s. As a consequence, the step-wise influencing path p (Definition [s]), which is the backbone of 
the compactification protocol, can not be defined anymore for these qubits. 

In summary, if one wants to save memory in the circuit model, the substitutions in Equations [T9] 



and|20]must be avoided, allowing the identification of all step- wise influencing paths and, consequently, 
the removal of auxiliary qubits. Conversely, if J-gate parallelization is the goal. Algorithm [3] can be 
easily adapted to create just \0^\ — p many J-blocks in the extended circuit (instead of \0^\ many), 
where p is the number of Pauli measurements in the associated measurement pattern. 



Consider the circuit shown in Figure 17 -a. It has two J gates with arbitrary angles (arbitrary 
angles are omitted in the Figure) and two Pauli angles, namely and |. Using the method in [TT] , 
we can translate this circuit to the MBQC model; the associated graph is shown in Figure [T7]-b and 
the flow measurement pattern is given by: 

Z'q^ X'q' X^' MJ/^ Z^^ Mf 1 Eg (21) 

where. Eg = E^qE4^E25E24E23Ei2. Using the signal shifting technique explored in this paper, we 
obtain the following SSF measurement pattern: 

Now let us see what happens when we use the fact that some measurements are Pauli measurements. 
Using Equation [T9j we can optimise it by removing the dependency between measurements 1 and 2, 
obtaining the following measurement pattern: 

Zg"^+"^Z3"^Xg"^+"^X3^2+"^M5°X5^^+^^M2''/^M4^^Mf^£;c (23) 

Alternatively, the measurement on qubit 5 can also be parallelised with the ones in 1 and 4, using 
Equation [20| The optimized measurement pattern is the following: 
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Figure 17: Using compactification procedures to obtain J-gate parallelization for circuits with Pauli 
angles. See example 2 for more information. 



Finally, both optimisations can be considered together. In this case, all measurements can be 
performed at once: 

Z^'^'^Z^'XI'^'^XI'^'^^'^M^^^M^M^^M^'Eg (25) 

With few adaptation^ Algorithm [s] can be applied to the extended circuits associated to the 
measurement patterns in equations [22] to [25| The obtained circuits are shown in Figures [TTl -c to 
17-f, respectively. Note that the optimization given by equations 19 and 20 in the MBQC model 
becomes J-gate parallelization in the circuit model. However, it is not clear in which cases the 
J-gate parallelization implies time optimization. An in-depth analysis of the time-memory tradeoff 
for extended circuits with pauli angles remains as an interesting open question. 



4 Complexity analysis 

In this section we analyse the space, i.e. the number of wires, and depth complexity of our proposed 
optimization scheme. We show that our procedure can give a more optimal circuit in terms of both 
space and depth compared with the general method obtained in [IJLJ. An overview of this analysis is 
shown in Figure [TSj 

Let C be our initial quantum circuit with n wires, m = poly(n) J gates and depth dc = poly(n). 
We can translate the computation implemented by C to an MBQC pattern P with m + n qubits, out 
of which m are measured qubits, and depth dp < d fTT]. This pattern can be further optimised to 
depth dss ^ dp by performing signal shifting [llj and obtaining a new pattern Pss- The extended 
quantum circuit Css, corresponding to the signal shifted pattern Pss can be created via the method 
given in Definition [TOj This new circuit Css performs exactly the same computation as the initial 
circuit C and has the same number of J-gates, m -\- n wires and depth dss • 0{m). One could apply 



^Due to the modifications introduced by equations 



19 



and 



20 it is not possible to create J-blocks for p many wires, 



where p is the number of times equations 19 or 20 were used in the original measurement pattern, 
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the parallelisation method of [10] to create a parahehsed circuit Cpar with depth dss • O(logn) and 
size 0{m?) [llj. Depending on the value of (iss^ the depth of Cpar can be smaher than that of the 
original circuit C, but this could increase the space used, which can increase considerably, i.e from 
n to 0{im?). 
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Figure 18: A summary of the optimisation procedure using extended translation and compacti- 
fication. The parallelisation procedure described in [llj (without Pauli simplification) ends with 
the "Parallelised Circuit" whereas our optimisation can end with a "Compactified Circuit" or a 
"Parallelised Compactified Circuit" depending the goal of the optimization. 

To resolve this problem, instead of parallelising the circuit Css we apply our compactification 
method. This will give us the compact circuit Ccom with n wires and dss • 0{deg{G)) where deg{G) 
is the degree of the open graph (G, /, O) of Pss- 

Lemma 28. Let Css be the quantum circuit obtained through extended translation from a signal 
shifted measurement pattern Pss on an open graph (G, /, O) with depth dss- The application of the 
compactification Algorithm^ to Css results in a circuit with depth dss ' 0{deg{G)). 



Proof. Algorithm [s] is the application of Rewrite Procedures (1) - (5) as defined in Section 3.2 to the 
extended circuit of Pss- Note that the number of J-gate layers in Css is equal to the depth of Pss 
[Hj. Since the rewrite procedures do not change the J-gate layers, the compactified circuit Ccom 
will also have dss J-gate layers. 

Now we will count the number of gates in-between the J-gates. We do not need to be concerned 
with gates in the C layers, since Algorithm [3] removes all of them. It is possible that the rewrite 
procedures leave some Eij — CZij gates in front of qubit i. As all of these gates correspond to an 
edge in the graph G, the maximum number of such gates has to be 0{deg(G)). According to the 
rewrite procedures, none of the CXj^^i gates created will be moved past the J-gates. As can be seen 
from the rewrite procedures and Algorithm |3) these will be created because of the existence of some 
T^f-^{i)k gates in Css- Since these correspond to the edges connected to the vertex f~^{i) there can 
only exist 0{deg{G)) many CX gates between any two J-gates. Now we know that there can be a 
total of 0{deg(G)) two qubit gates in-between any two J-gates and hence also between the J-gate 
layers. Therefore the total depth of the circuit Gcom will be dss ' 0(deg(G)) □ 
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Already we can see, that if deg(G) < logn the Ccom wih have smaller depth than Css while 
the corresponding space will be considerably smaller. We can further decreased the depth of Ccom 
by applying the parallelisation method from [lOj. The depth of the new circuit CcomPar will be 
dsS'0{log{degG)) and size O(n^). Note that because of the way C is translated into P, the maximum 
number of edges connected to a vertex in G will be 0((i), where d = poly(n). Therefore the depth of 
CcomPar Can be written as dss ' O(logn). Hence the compactification procedure together with the 
parallelisation method from [lOj will in the worst case give us the same depth as the method from 
pT] , without Pauli Simplifications, but uses considerably less qubits {n? vs m^, where m = poly{n)). 

5 Discussion and summary of results 

Initially MBQC was proposed as an alternative architecture for the implementation of quantum 
computing. However from early on the distinct parallel power of the model attracted researchers 
to explore further this unique feature of the model. In a series of results the key concepts of flow, 
signal shifted flow, gflow, maximally delayed gflow, focused gflow and information preserving flow 
were introduced flSl IH [3 El [29j. They address the general question of determinism in MBQC 
while shedding light on the parallelism as well. Although it is now known that the MBQC parallel 
power is equivalent to the quantum circuit with unbounded fanout [30j, further investigation is 
required to fully take advantage of this extra power. In this paper we continue this line of research 
by first presenting a surprising link between signal shifted flow and maximally delayed gflow. The 
surprise comes from the fact that the former is obtained via a simple rewrite rules of pushing the Z 
dependencies of a pattern to the end of the computation, while the latter is constructed directly 
from the stabilisers of the underlying graph. This leads to a new efficient procedure for finding the 
optimal gflow of graphs with flow as we discussed in the last section. 

Moreover the link between signal shifted flow and optimal gflow opens a new direction to unify 
further the "flow" structure and fully characterise the constructions behind the parallel power of 
MBQC. To begin with, we succeeded in extending the applicability of the compactification method 
from [13j by proving that it can also be used for the signal shifted flow. This will allow us to translate 
back the obtained parallel structure of the MBQC into the circuit model without increasing the 
size, leading to an automated optimisation procedure for quantum circuits. The automated scheme 
explore the global structure of circuit to parallelise several J-gates and to group together several 
CX gates, which allow further optimisations to take place using specific method for Clifford gate 
parallelisation. Interestingly the scheme fails to compact the parallel pattern obtained via Pauli 
simplification rules [11]. In other words one needs to keep the extra space to keep the parallel 
depth obtained due to the Pauli measurements. This further indicates the crucial role that Clifford 
computation (corresponding to Pauli measurements) play in obtaining the superior parallel power of 
MBQC over quantum circuit. 

The above trade-off naturally suggests the consideration of a hybrid model for quantum computing, 
where part of the computation is processed using the quantum circuit model and the other part 
using MBQC. Such a computation can be obtained by a partial compactification of the extended 
circuits. First a quantum circuit acting on n qubits is translated to the MBQC model, where it 
requires m > n qubits. Then, we use an automated compactification procedure to obtain many 
different circuits implementing the same computation but with different depth and number of qubits 
(ranging from n to m many qubits). In those circuits, part of the computation is implemented using 
MBQC to obtain the parallel depth due to Pauli measurement while the rest are performed in the 
circuit model as one would not achieve any parallel advantage. We believe that this can be of great 
value for experimental implementation of quantum computation, since one can design algorithms 
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more adapted to the available experimental resources. 

The application of the different techniques introduced in this paper to known quantum algorithms, 
as well as a full comparison with other known optimization methods beyond MBQC, constitute an 
interesting subject to be explored in the future. 

Acknowledgments 

We are grateful to Ernesto F. Galvao for many helpful discussions. We acknowledge financial 
support by the Instituto Nacional de Ciencia e Tecnologia de Informacao Quantica (INCT- IQ/CNPq 
- Brazil) that sponsored R. Dias da Silva visit to Scotland where this work was completed. E. 
Kashefi also acknowledges support from UK Engineering and Physical Sciences Research Council 
(EP/E059600/1). 

References 

[1] D. Deutsch. Procedmgs of the Royal Society 425 73-90 (1989). 

[2] M. A. Nielsen and I. L. Chuang. Quantum Computation and Quantum Information (Cambridge 
University Press, 2002). 

[3] R. Raussendorf and H. J. Briegel. Phys. Rev. Lett. 86, 5188-5191 (2001). 

[4] V. Danos, E. Kashefi and P. Panangaden. Journal of the ACM 54, 2 (2007). 

[5] D. E. Browne, E. Kashefi, M. Mhalla, and S. Perdrix. New J. Phys. 9, 250 (2007). 

[6] R. Raussendorf, H. J. Briegel, and D. E. Browne. Journal of Modern Optics 49, 1299 (2002). 

[7] J. C. Garcia-Escartin and P.Chamorro-Posada. arXiv:1110.2998vl [quant-ph] (2011). 

[8] M. Sedlak and M. Plesch. Central European Journal of Physics, 6, 128-134 (2008). 

[9] D. Maslov, G. W. Dueck, D. M. Miller and C. Negrevergne. IEEE Transactions on Computer- 
Aided Design of Integrated Circuits and Systems 27(3), 436-444 (2008). 

[10] C. Moore and M. Nilsson. SIAM J. Computing. 31, 799-815 (2001). 

[11] A. Broadbent and E. Kashefi. Theoretical Computer Science 410 (26), 2489 (2009). 

[12] R. Duncan and S. Perdrix. Proc. of ICALP Part II, 6199, 285-296 (2010). 

[13] R. Dias da Silva and E. Galvao. arXiv: 1209.5079 [quant-ph] (2012) (2012) 

[14] M. Mhalla and S. Perdrix. Proc. of 35th ICALP, 10, 857-868 (2008) 

[15] V. Danos and E. Kashefi. Phys. Rev. A 74, 052310 (2006). 

[16] D. Gottesman and I. Chuang. Nature 402, 390-393 (1999) 

[17] S. F. Huelga, M. B. Plenio and J. A. Vaccaro. Phys. Rev. A, 65, 042316 (2002). 

[18] S. F. Huelga, J. A. Vaccaro, A. Chefles and M. B. Plenio Phys. Rev. A, 63, 042303 (2001). 

[19] D. Leung. arXiv:0111.122 [quant-ph] (2001). 



49 



[20] M. Nielsen. Phys. Lett. A. 308 (2-3): 96-100 (2003) 
[21] R. Jozsa. arXw:0508.124 [quant-ph] (2005) 

[22] R. Raussendorf, D. E. Browne and H. J. Briegel. Phys. Rev. A 68, 022312 (2003). 

[23] M. Hein, W. Diir, J. Eisert, R. Raussendorf, M. Van den Nest, and H. J. Briegel. In Proceedings 
of the International School of Physics "Enrico Fermi" on "Quantum Computers, Algorithms and 
Chaos" 162 (2006). 

[24] D. Gross, S.T. Flammia and J. Eisert. Phys. Rev. Lett. 102, 190501 (2009). 

[25] M. Van den Nest, A. Miyake, W. Diir, and H. J. Briegel. Phys. Rev. Lett. 97, 1505 (2006). 

[26] C. E. Mora, M. Piani, A. Miyake, M. Van den Nest, W. Diir, and H. J. Briegel. Phys. Rev. A 
81, 042315 (2010). 

[27] T.-C. Wei, I. Affleck, and R. Raussendorf. Phys. Rev. Lett. 106, 070501 (2011). 
[28] R. Duncan. arXiv: 1203.6242 [quant-ph], (2012). 

[29] M. Mhalla, M. Murao, S. Perdrix, M. Someya and P. S. Turner. TQC (2011). 

[30] D. E. Browne, E. Kashefi and S. Perdrix. Lecture Notes in Computer Science 6519, (2011). 

[31] N. de Beaudrap. Phys. Rev. A 77, 022328 (2008). 



50 



